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The theoretical effort includes the design of sampled 
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I. INTRODUCTION 



Signal processing covers a wide range of application in 
such diverse fields as biomedical engineering, acoustics, 
sonar, radar, speech communication, data communication, and 
many others. As an example, one may wish to remove interfer- 
ence, such as noise, from the signal or to modify the signal 
to present it in a form which is more easily interpreted by 
an expert. As another example, a signal transmitted over a 
communications channel is generally perturbed in a variety of 
ways, including channel distortion, fading, and the insertion 
of background noise. One of the objectives at the receiver 
is to compensate for these perturbances. In each case, 
processing of the signal is required. 

Signal proce-ssing was carried out typically by using 
analog equipment before the advent of digital computers which 
promoted the intensive and successful development of digital 
signal processing in the past decade. At present, much of 
the digital signal processing is performed on general purpose 
digital computers using software. However, tremendous progress 
has also been made in other types of signal processors using 
electronic firmware and hardware. Particularly, using elec- 
tronic technologies, a new class of sampled analog processors 
is being developed using charge transport devices. Their 
signals are analog but the independent variables are discrete 
which can be either time variable, spatial variables or other 
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transform variables. The signal processing implementations 
are illustrated in Table I. I. 

Signal processors can be grouped into two classes: 
filtering and spectral analysis. This classification is 
highlighted in Figure 1.1. All of these signal processing 
systems require the basic essential functions, such as cor- 
relation, convolution and transformation (Fast Fourier trans- 
form). In other words, they are shown to be centered around 
four basic mathematical operations: delay (or shift), mul- 

tiplication, summation and generation of special functions. 
These operations have been performed by digital techniques 
and it was probably the most convenient method. However, 
the new sampled analog signal processing systems only require 
a storage for discrete time analog samples and not discrete 
or quantized signal amplitudes; in fact, because of the 
computation time, frequency selectivity, and part cbunts, it 
is undesirable to quantize the signal. Therefore, discrete 
time analog devices are being developed to perform the convo- 
lution, correlation, etc. , without the need to convert the 
signal to digital format. 

There are several distinct advantages of using the discrete 
time analog devices over the digital techniques to perform the 
same algorithm: 

(1) In general, the discrete time device system can perform 
the same computation with much greater speed over the 
digital technique. 

(2) The device stores discrete time analog amplitudes; 
therefore, it does not have the quantization errors 
inherent in the digital system. 
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(3) The actual multiplying and summing operations are 
easily implemented, and with fewer components than 
the digital system. 

To illustrate the advantages that these discrete time 
devices have over the digital technique, a simple system 
shown in Figure 1-2 is used. This is a model of a simple 
single pole, first order discrete time filter. To implement 
this system with a discrete time device such as the SAD (Serial 
Analog Delay) , all that is required is a straight substitution 
of the analog device for the delay block, an operational 
amplifier for the adder and a potentiometer proportional to 
the weighting coefficient. This system is shown in Figure 
1-3. 

Comparing this system implementation to the digital form 
in Figure 1.4 demonstrates the difference between the two 
systems. It requires at least two memory registers, one to 
store the sampled data, the other to store the multiplying 
constant. In the figure it is seen that more devices are 
required, such as the adder and the costly analog to digital 
converters, as well as the digital to analog converter to 
complete this system. 

The integrated circuit delay devices can perform in 
principle most of the necessary algorithms used in most equip- 
ment for signal processing today. However, due to the primi- 
tive nature of these new classes of devices, there are 
limitations in some fields of application. 

The purpose of this study is to perform the evaluation of 
Reticon TAD-12 tapped delay lines available commercially and 
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to find out their performances and limitations in order to 
use them properly. 



TABLE I . I Signal Processing Implementation 



Signal 



Equipment 



Implementation 



Analog 



Digital 



Sampled 

Analog 




General Purpose Computer 
Special Purpose Computer 
Mini/Microcoraputer 
(Microprocessor) 

Signal Processing Function 

Modules : 

Multiplier 

Correlator 

Convolver (matched 
filter) 

Filter 

Discrete Fourier 
Transformer 

Frequency synthesizer 

etc. 



Software 

Firmware 



Hardware 
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FIGURE 1.1 Classification of Signal Processors 
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FIGURE 1 . 3 Equivalent Circuit Implementation 




FIGURE 1 .4 Digital Equivalent Circuit Implementation 
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II. SAMPLED ANALOG NON- RECURSIVE FILTER 



A. PRINCIPLE 

The transfer function of a digital or sampled analog 
filter can be expressed, in Z domain, by one of the follow- 
ing equations. 



Recursive filter: H(Z) = 

1 * h„ z‘" 

1 m 

Non-Recursive filter: H(Z) “ Zq ^ ^ 



C2-1) 

C2-2) 



The block diagram of a sampled analog recursive filter 
is illustrated in Figure 2-1: 



E. (z) 



m 







FIGURE 2.1 Schematics of a Recursive Filter 
(Tapped Delay Line) 
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E 



in 




E + 

out 

FIGURE 2,2 Schematic of a Non-Recursive Filter 
(Tapped Delay Line) 



In the case of causal system with finite duration of 
impulse response, the system function can be written in the 
form: 



H(Z) 



N--1 




_n 

h(n) Z 



(2-3) 



{h(n)} is the impulse response defined over the time interval 
0 £ nT £ (N - 1)T, H(Z) is a polynomial in Z ^ of degree N - 1. 
Thus, H(Z) has N - 1 zeroes in the finite Z- plane. 

The frequency response H(e^^) is obtained by substituting 
Z by e^^^. 
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(2-4) 



jwT. N-1 -jwrtT 

H(e*' ) = Z h(n) e 

n=0 

We recall that any finite duration sequence is completely 
specified by N samples of its Fourier Transform [1] , so that 
the design of a Finite Impulse Response Filter may be 
accomplished by finding either its impulse response coeffi- 
cients or N samples of its frequency response. Both methods 
are discussed in the following section. 



B. LINEAR PHASE FINITE IMPULSE RESPONSE FILTER 

In many applications, e.g., speech processing, data 
transmission, it is desirable to design filters to have 
linear phase. In this way, signals in the passband of the 
filter are reproduced exactly at the filter output except 
for a delay corresponding to the slope of the phase. One 
of the most important features of FIR systems is that they 
can be designed to have exactly linear phase. The unit 
impulse response for a causal FIR system with linear phase 
has the property that: 

h(n) = h(N - 1 - n) (2-5) 



To see that this condition implies linear phase, we write 
equation 2-4 as: 



H(2) = 



(N/2)-l 

Z 

n=0 



(N/2)-l 

Z„ 



h(n) Z'^ + 
h(n) + 



N-1 

-z- 

n=N/ 2 



(N/2)-l 

Z 

n=0 



h(n) Z 

h(N-l-n) z" (N-1-^) 



25 



where N is assumed to be even. Using equation 2-5, one can 
write : 

-fN-l-n'l 

H(Z) = I h(n)[Z n + 2 (2-6) 

n=0 



If N is odd, it is shown that 



Evaluating equations 2-6 and 2-7 for Z = e^'^, one obtains, 
for N even, 

H(e>) = ,-jw[(N-l)/2](W/’-l 2 h(„) cos [w(n - 

n=0 ^ 

(2-8) 

and for N odd. 






[(N-3)/2]-l , 

+ Z 2h(n) cos [W(n - ^V^) ] } (2-9) 

n=0 ^ 



In both cases the sums in brackets are real, implying a 
linear phase shift corresponding to delay of (N-l)/2 samples. 
For N odd, the phase shift corresponds to an integer number 
of samples delay. For even N, the delay is an integer plus 
one-half of sampling period. This distinction between odd 
and even values of N is often of considerable importance in 
design and realization of FIR filter. Examples of impulse 
response having linear phase are shown in Figure 2-3. 
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FIGURE 2.3 Typical Impulse Response for Linear Phase FIR Filters 

C. DESIGN TECHNIQUE OF FIR FILTERS 

1 . Determination of Impulse Response 

Since the frequency response of H(e^^) of any digital 
or sampled analog filter is periodic in frequency, it can be 
expanded in Fourier series of the form: 

HCe^'^) = h(n) e’^^” (2-10) 

n=-“ 

where the sequence h(n) playing the role of the ’’Fourier 
coefficients,” i.e., 

h(n) = ^ HCe^"") dw (2-11) 



27 



One can easily recognize that the coefficients of the 
Fourier series, h(n) , are identical to the impulse response 
of a digital filter. There are two difficulties with the 
representation of (2-10) for designing FIR filters. First, 
the filter impulse response is infinite in duration since the 
summation in (2-10) extends to ±o° . Second, the filter is 
unrealizable because the impulse response begins at -«>; i.e., 
no finite amount of delay can make the impulse response 
realizable . 

One possible way of obtaining an FIR filter that 
approximates H(e^^) would be to truncate the infinite Fourier 
series (2-10), at n = ± M. But direct truncation of the 

series leads to the well known Gibbs phenomenon. However, 
the non-uniform convergence phenomenon can be moderated 
through the use of a less abrupt truncation of the Fourier 
series by applying the so-called "windowing” or "weighting 
function. " 

2 . Effect of Window 

A more successful way of obtaining an FIR filter is 
to use a finite weighting sequence W(n) , called window, to 
modify the Fourier coefficients h(n) in (2-10) to control the 
convergence of the Fourier series. The technique of windowing 
is illustrated in Figure 2-4. At the top of this figure is 
shown the desired periodic frequency response H(e'^'^) and its 
Fourier series coefficients (h(n)}. The next row shows a 
finite duration weighting sequence W(n) with Fourier trans- 
form W(e^^) . W(e^^^) , for most reasonable windows, consists 
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of a central lobe which contains most of the energy of the 
windows and side lobes which generally decay rapidly. To 
produce an FIR approximation to H(e^^), the sequence 



r. _ j.h(n) . W(n) 
h(n) = {q 



-M < n < M 
elsewhere 



( 2 - 12 ) 



is formed. The third row of Figure 2.4 shows h(n) and its 
Fourier transform H(e^^) , which is the convolution of 
E:(e^^) and W(e^^) , since h(n) is the product of the sequences 
h(n) and w(n) . The last row of Figure 2-4 shows the realiz- 
able sequence g(n) , which is a shifted version of h(n) and 
may be used as the desired filter impulse response. 

As seen in the example of Figure 2-4, there are 
several effects on the resulting frequency response of win- 
dowing the Fourier coefficients of the filter. A major effect 
is that discontinuities in HCe^'^) become a transition band 
between values on either side of the discontinuity. Since 
the final frequency response of the filter is the circular 
convolution of the ideal frequency response with the window's 
frequency response, it is clear that the width of these tran- 
sition bands depend on the width of the main lobe of W(e^'^) . 

A secondary effect of windowing is that ripple from the side 
lobes of WCe^'^) produces approximation errors (ripple in the 
resulting frequency response) for all w. 

There are many windows proposed that approximate the 
desired characteristics. But, in general, the desirable 
window should have the following characteristics: 
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FIGURE 2.4 Illustration of windowing 
Desired periodic frequency response 
h(n) Its Fourier series coefficients 
w(n) Finite duration weighting sequence 
W(e^^) Its Fourier transform 
h(n) = h(n) . w(n) 

H(e^’). » W(e^’) 
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1. Small width of main lobe of the frequency response of 
the window containing as much of the total energy as 
possible . 

2. Side lobes of the frequency response that decrease in 
energy rapidly as w tends to it. 

In the design of the filters in this study, the 

Hamming window is chosen. It is specified by the equation: 

, 0.54 - 0.46 cos 0 < n < N-1 

w(n) = { — — 

0 , elsewhere (2-13) 

Figure 2-5 illustrates the Hamming window and the rectangular 
window. 



W(n) Rectan^lar 




FIGURE 2.5 Hamming Window and Rectangular Window 

3 . Frequency Sampling Technique 

Figure 2-6 shows an arbitrary frequency response 
(solid curve) that one wants to approximate and a sequence of 
N frequency samples H^(Figure 2-7) that can be represented as 

H(k) = |H(k)|ej®^^\ k = 0, 1..., N-1 
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N 

FIGIffiE 2.6 Desired Continous Frequency Response 




FIGURE 2.7 Frequency Samples 



Using the inverse discrete Fourier transform, a 
finite duration impulse response can be determined from H(k) 



h(n) 



H(k) eH2’'/N)nk 

^ k=0 



(2-14) 



where H(k) can also be written in the form: 



H(k) = h(n) 
n=0 



(2-15) 



D. APPLICATIONS 



1 . Bandpass Filter Design Algorithm 

In the analog or continuous time system, the bandpass 
filter has the transfer function H(f). 




FIGURE 2.8 Bandpass filter frequency response 

A, f, < f < f. 

H(f) = { ^ ^ 

0, otherwise 

In the discrete time system, a rectangular bandpass 
filter has the frequency response H(e^^) as shown in Figure 
2-9. It is periodic of period 2 it . 




0, otherwise. 



Ws 

with €2 - sampling frequency. 
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A 




PIGDEE 2.9 Discrete Time Bandpass Filter Frequency Response 



Basing upon the periodic nature of a discrete time 
filter frequency response, the impulse response can be com- 
puted by finding the Fourier coefficient, using equation 2-11. 

h(„) . f A eJ"" dw (2-16) 

^TT VT 

1 

1 ^2 

= ± f A (cos wn + j sin wn) dw. 
ir w 

1 

Since only real term may be implemented, the series coeffi- 
cients are: 

Case I : n ^ 0 

1 ^2 

h(n) = i A cos (wn) dw 

= — (sin Zirf - sin Zirf ) (2-17) 

irn ^ 2 1 

Case 2: n = 0 

Let w = 2 it (fo + 

W = 2 7T (fo + 

One can write using equation (2-16): 
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ur ^ A .^ 2 - jwn J 

h(n) = — / dw 

ir Wj 

A r jw^n jw,riT 
" “ - eJ . ] 

” t4t (sin ;r n4 f) 

j nTT ^ 

= 2A A f (cos 2Trnfo + j sin 2 tt n £o) 

(2-18) 

When n tends to 0, the real part o£ equation (2-18) tends 
to: 

h(0) = 2AA£(1.)(1.) = 2AA £ (2-19) 

Window weighting £unctions are used to modi£y the 
coe££icients in order to improve the £ilter characteristics. 

I£ the Hamming window is used, the weighted coe££icients 
A(n) become 

A(n) = h(n)(0.54 + 0.46 cos (^tt)) (2-20) 

j 

A computer program shown in Appendix C was established 
to design the non-recursive bandpass £ilter. It was based 
upon the equations 2-17, 2-19, 2-20. 

2 . Prewhitening and Dewhitening Filters in 

Communication System 

The prewhitening and dewhitening £ilters are also 
called preemphasis and deemphasis £ilters in the communication 
£ield. 

Suppose one desires to transmit a baseband signal 
using FM modulation and requires the best possible signal to 
noise ratio. One possible solution is to raise the level o£ 
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modulating baseband signal to the maximum extent possible 
in order to modulate the carrier as vigorously as possible. 

But the modulating signal level may be raised only until the 
distortion exceeds a specified value (maximum allowable value) . 

However, the baseband signal happens to be an audio 
signal; it turns out that something further can be done. An 
audio signal usually has the characteristic that its power 
spectral density is relatively high in the low frequency 
range and falls off rapidly at higher frequencies. For ex- 
ample, speech has little power spectral density above about 
3 kHz. Music extends farther into the high frequency range, 
the feature still persists that most of its power is in the 
low frequency region. As a consequence, when one examines 
the spectrum of the sidebands associated with a carrier which 
is frequency modulated by an audio signal, one finds that the 
power spectral density of the sidebands is greatest near the 
carrier and relatively small near the limits of the allowable 
frequency band allocated to the transmission. The manner in 
which one may take advantage of these spectral features, 
which are characteristic of audio signals, in order to improve 
the performance of an FM system is shown in Figure 2-10. 

At the transmitting end, in Figure 2-10, the baseband 
signal m(t) is not applied directly to the FM modulator but 
is first passed through a filter of transfer characteristic 
Hp (w) , so that the modulating signal is mp(t). The modifi- 
cation introduced into the baseband signal by the first 
filter is undone by the receiver filter which follows the 
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discriminator and has transfer characteristic 1/Hp(w) . The 
noise passes through only the receiver filter which may be 
used to suppress the noise to some extent. 

The selection of the transfer characteristic Hp(w) 
is based on the following considerations. At the output of 
the demodulator, the spectral density of the noise increases 
with the square of the frequency, as shown in Figure 2-11. 

It is given by the equation: 

Gn(f) “ ^ ■ |f I 1 I (2-21) 

where a = constant proportional to limited amplitude 

of the carrier obtained with the hard limiter 
in the demodulator 
A = amplitude of the carrier 

n = power spectral density of the white noise 

B 

over the range |f| < y • 




FIGURE 2.11 Power Spectral Density at the Output 
of an FM Demodulator 
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Hence the receiver filter will be most effective in 
suppressing noise if the response of the filter falls off 
with increasing frequency. In such a case the transmitter 
filter must exhibit a rising response with increasing fre- 
quency. The premodulation filtering in the transmitter, to 
raise the power spectral density of the baseband signal in 
its upper frequency range, is called preemphasis. The fil- 
tering at the receiver to undo the signal preemphasis and 
to suppress noise is called deemphasis. 

Referring to Figure 2-10, one requires that the nor- 
malized power of the baseband signal m(t) must be the same 
as the normalized power of the preemphasized signal m^Ct). 

If Gjjj(f) is the power spectral density of m(t) , the density 
of nip(t) is |Hp(f)l^ Gjjj(f) , and one requires that: 

fm fm 

where 

fm = maximum frequency of modulating signal. 



In the absence of deemphasis, output noise is: 



No = ^ 

With deemphasis, the output noise is: 

fm 



f'od ■ “’’'I linTT 



df 



(2-23) 
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IS 



The ratio, R 




> 



R 



fmV5 



(2-24) 



Since the signal itself is unaffected in the overall 
process, the quantity R is the ratio by which preemphasis- 
deemphasis improves the signal to noise ratio. 

The prewhitening filter has the transfer character- 
istic illustrated in Figure 2-12. The dewhitening filter 
has the transfer characteristic illustrated in Figure 2-13. 
The relation between these two filters is: 

HdCs^") = (2-25) 

“ Hp(eJ") 

where 

= transfer function of prewhitening filter 
Hp(e^^) = transfer function of dewhitening filter. 

3 . Prewhitening Filter Design Algorithm 

a. In the analog or continuous time system, the 
prewhitening filter has the transfer characteristic shown 
in Figure 2-12 and has the transfer function. 



H (f) 
P 



R , |£1 < fo 

R * A(f - £o) , fo i |£| 1 fp 

(2-26) 
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FIGURE 2.13 



f- f f nonnalized 

0 p 

Dewhitening Filter Transfer Function 



b. In the sampled analog or digital system, the 
frequency response of the prewhitening filter is periodic. 
It is illustrated in Figure 2-14. 
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Analog Prewhitening Filter 



The frequency F is normalized with respect to 
the sampling frequency . 

The design of a linear phase prewhitening filter 
consists of the calculation of the impulse response which 
can be computed using (2-11) 



h(n) = i dw 

i II 

which becomes, using ^ = F 



h(n) 

h(n) 



2 R dF + 2/p^ [R + A(F-Fq) ] e^ ^’^”^dF 

^o 



U 



V 



(2-27) 



Since h(n) is real, one can compute only the real 
part of U and V. 
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Case 1: n ^ 0 



> 



F = ^ 



R 

Re[U] = 2R cos (2TrnF)dF = — sin 2irn Fq 



o 

F, 



Re[V] = 2R / ^ cos (2TrnF)dF + 2A / ^ (F - Fq) cos (2TrnF)dF 






+ 2A (cos 2iTnFp - cos 2iTnFo) 

F_ Fr) 



AFq 

irn 



[sin 2irnFo - sin 2irnFp] 



h(n) = Re [U] + Re [V] = 2irn Fq ] 

(2-28) 

Case 2 : n = 0 , ^p ~ ^s/2 



With n = 0, (2-27) becomes: 

Fp Fp 

h(o) = 2/ RdF + 2/ [R + A(F - Fq) ]dF 
op 
^0 

hCo) = [R * (Fp - Fo)=l 



2Fp 



2Fp '-"P 



(2-29) 



Taking into account the constraint of the linear 
phase FIR filter: 

h(n) = h (N - 1 - n) 

and using (2-28) and (2-29), a computer program 
is prepared to calculate the impulse response of 
the filter, a Hamming window function is used to 
improve the amplitude and phase response of the 
prewhitening filter. 
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4 . Dewhitening Filter Design Algorithm 

In the design of the complementary dewhitening filter, 
one can use the frequency sampling technique. Using (2-15) 
and the given impulse response of the prewhitening filter, 
the frequency samples 

{Hp(k) } , k = 0, 1 , N-1 

can be computed and stored in the memory of the computer. 

Then the frequency samples of dewhitening filter 
Hjj(k) can be computed as: 

^D ~ H Ck)' * k=0, 1,....,N-1 

Using the inverse discrete Fourier transform given 
by (2-14) , the finite impulse response of the dewhitening 
filter can be obtained. 

Using an appropriate window, such as Hamming window, 
the frequency response of the dewhitening filter can be 
improved also. 
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III. EXPERIMENTAL EVALUATION OF TAD- 12 



A. SYSTEM CONFIGURATION 

Reticon Corporation markets a Tapped Analog Delay Device 
which provides 12 output taps with the following time delay 
configuration: 




* PIGDHE 5*1 ' BqvilTalent Circuit of TAD- 12 

The first delay section is delayed by one clock period, 
the rest have delay time of two clock periods. 

The delay lines are structurally organized as shown in 
Figure 3-2. One delay line contains N capacitive storage 
elements with two I/O access switches. One access switch 
connects the storage element to a common input line and the 
second switch connects it to the common output line. One 
register and its associated multiplexing switch sequentially 
time samples the analog input signal onto the ”N" capacitive 
storage elements. Information is read onto the nth element 

1“ H 

while simultaneously reading out the n^ 



+ 1 element. 



CLOCK INPUT 







FIGURE 5»2 Organization of a Single Delay Line 



46 



Each element has an output buffer amplifier to provide 
a zero order hold of the sampled signal. This produces a 
sampled and held output. 

The test circuit for TAD-12 is shown in Figure 3-3. It 
consists of a pulse shaper circuitry and a DMOS flip-flop. 
Because the tap outputs are at +5 volts bias level, due to 
the TAD-12 circuitry, the resistor 2 kilo-ohms of the oper- 
ational amplifier non- inverting input is connected to a +5 
volts rather than ground. The output of the summer, there- 
fore, rides on a 5 volts bias level. 

B. FREQUENCY DOMAIN EVALUATION 

Using the circuit board shown in Figure 3-3, the frequency 
dependance, the non-uniformity, the harmonic distortion are 
evaluated at different sampling frequencies for three differ- 
ent tap resistances: 

= 10 Kfi , 100 Kfi , 1 Mfi . 

1 . Frequency Roll-off 

The frequency characteristic of the sampled and hold 
process is of the ^ shape with x = ^ . 

X tg 

At sampling frequency f = 500 KHz, for R = 10 Kfi, 

^ L. 

the roll-off is close to the theory. When the tapping 
resistances R^.^p is increased, the roll-off is more serious 
and wriggles around. 

Figure 3-4 illustrates the magnitude of the frequency 
response of Tap #1, Tap #6 and Tap #11 outputs using 10 kilo- 
ohms resistor. The solid curve represents the theoretical 



response . 
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FIGURE 5* 5 Test Circuit for Reticon TAD -12 Evaluation 
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+5V -15V +15V 



From this experiment it is found that the TAD- 12 

output has rather complicated frequency variations although 

s in X 

they follow the theoretical relation roughly if the 

tapping resistor is 10 Kilo- ohms and the sampling frequency 
less than 2.5 MHz. The loading effect is much more pronounced 
if the tapping resistors are equal or greater than 100 Kilo- 
ohms . 

2 . Non-Uniformity (Spatial or Fixed Pattern Noise) 

The non-uniformity of different taps outputs varies 
depending on the tap resistance, Rj^, the sampling frequency, 
f s , and the input signal frequency, f. The general result 
at sampling frequency of 500 KHz is the following: 

- For f < 10 KHz, the non-uniformity is within 10% 

- For f > 10 KHz, the non-uniformity can deteriorate 
to 15% - 50% range. 

Figures 3-10 through 3-17 illustrate the non-uniformity 
result of sampling frequency of 500 KHz, for different input 
signal frequencies and for different tap resistors: 

= 1 Kj^, 10 Kn, 100 Kn, 1 Mn. 

At high sampling frequencies, using tap resistances 
less than 100 Kn, the uniformity of tap output can be improved 
at higher input signal frequencies. 

At fg = 0.5 MHz, Rj^ = 10 KJ^ and input signal frequencies: 

- f < 50 KHz, the non-uniformity is within 10% 

- f > 50 KHz, the non-uniformity varies from 13% to 25%. 
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At £5 = 2.5 MHz, Rj^ = 100 Kfi and input signal frequencies: 

- £ < 80 KHz, the non- uniformity is within 101 

- £ > 80 KHz, the non-uniformity varies from 11% to 70% 

At £g = 2.5 MHz, = 1 Mfi, and input signal frequencies: 

- f < 30 KHz, the non-uniformity is within 20% 

- £ > 30 KHz, the non-uniformity is greater than 30%. 

Figures 3-18 through 3-28 illustrate the non-uniformity 
data at sampling frequency equal to 2.5 MHz, for different 
typical tap resistors and for different input signal fre- 
quencies . 

Figures 3-24 through 3-40 illustrate the pictorial pres- 
entation of Reticon TAD- 12 performance at different sampling 
f re*quencies , using different tap resistances. 

3 . Temporal Noise Behavior of Different Taps 

The temporal noise of different taps is studied at 
the sampling frequencies higher or equal to 2.5 MHz. It 
varies, depending on the sampling frequencies and the tap 
resistances used. However, it should be pointed out that 
the tap #1 output has little noise at all sampling frequen- 
cies and tapping resistors used. It is shown in Figures 3-24, 
3-25, 3-26, 3-28, 3-40. The general temporal noise behavior 
of other taps is shown in Table III. I and Figures 3-24, 3-25, 

3-27, 3-35, 3-38. 
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TABLE III. I Noisiest Tap Outputs 



K 

tap 

1 Meg-ohm 100 Kilo-ohms 10 Kilo-ohms 

2.5 MHz Tap #10 Tap #11 Tap #12 

3.5 MHz Tap #10 Tap #4 Tap #2 

The temporal noise of different taps is not serious 
at low sampling rate (less or equal to 1.25 MHz). Figure 
3-41 illustrates the output of tap #1 and #12 at sampling 
frequency at fg = 1.25 MHz using 1 Meg-ohm tapping resistors. 

4 . Harmonic Distortion 

The second harmonic distortion, data is measured with 
the input frequency 10 KHz for different input amplitudes at 
different sampling frequencies: 

fs = 500 KHz, 1.25 MHz, 2.5 MHz 

and using different tap resistors: 

Rtap = 10 Kfi, 100 1 MJ^. 

The curves for the quantity 20 log y fundam^^al" 
are plotted and shown in Figures 3-42, 3-43, 3-44. The general 
results of the above measurement are summarized in Table 3. II. 

The loading effect on the dynamic range of the TAD-12 
is very pronounced at high sampling frequency (2.5 MHz). In 
order to get the second harmonic amplitude less than 1% of 
the fundamental, one has to use less than .12 volt rms input 
voltage if using 1 Mn tap resistor and less than 2 volts rms 
input voltage if using 10 Kn tap resistor. 
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FIGURE 3.4 Amplitude Frequency Response of Reticon TAD-12 Tapped Delay Line. 



0 DB 




00 

I I 



C\4 



VO 



I 



I 



53 



FIGURE 3*5 Amplitude Frequency Response of Reticon TAD - 12 Tapped Delay Line 



CO 

CO 

0 

O 

o o 




54 



FIGURE 5»6 Frequency Response of Reticon TAD-12 Tapped Delay Line Showing 
Its Bandwith, Nonuniformities Among Taps and Loading Effect 
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FIGURE 3»7 Amplitude Frequency Response of Reticon TAU-12 Tapped Delay Line « 
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FIGURE 3.8 Amplitude Frequency Response of Reticon TAD-12 Tapped Delay Line. 
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FIGURE 5.9 Amplitude Frequency Response of Reticon TAD-12 Tapped Delay Line, 



Number > i 
of Taps' 



20% vsuriation 



3 ‘ 
2 ' 



.04 



.045 

Tap Output Amplitude 



.05 volts 



PIGDEE 3.10 



Non Uniformity Data, 500 KHz 

R. - 1 Kilo-ohm. f. .» 45 KHz 
tap input 




FIGURE 3*11 Non Uniformity Data, f =* 500 KHz 

samp 

R. » 1 Kilo-ohm. f. .= 60 KHz 
tap input 
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PIGDEE 5.12 Non Uniformity Data, f^„„ - 500 KHz 
"tap" ^inpat’ 




Tap Output Amplitude 

PIGDEE 3.13 Non Uniformity Data. f„^„ - 500 KHz 

ScUiip 

E. =10 Kilo-ohms, 100 KHz 

tap input 
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Nimber » 
of Taps 



225^ variation 



5 • 

2 . 

1 .. 



.075 .1 •105 .11 .12 

Tap Output Amplitude 

FIGURE 3.14 Non Ifciformity Data. f„^ “ 500 KHz 

Scu&p 

“tap- Itilo-olms. 50 KHa 



Number 
of Taps 




FIGURE 3.15 Non Uniformity Data. 

R. » 100 Kilo-ohms, 
tap 



f = 
samp 

^input 



500 KHz 
•* 60 KHz 
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Number <■ 
of Tape 



11,59^ variation 



5 * 
2 • 



< » 






1 • 



.12 



.15 



i I 




Tap Output Amplitude 

FIGURE 5.16 Non l&iiformity Data, f = 500 KHz 

samp 

®tap" ■' ^input* ^ 



volts 




FIGURE 3*17 Non Uniformity Data, f = 500 KHz 

sanip 

R. » 1 Meg-ohm. f. .*120 KHz 
tap ° input 
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Tap Output Amplitude 

FIGURE 5 •18 Non Uniformity Data, f * 2.5 MHz 

samp 

R. - 10 Kilo-ohms. f. , 50 KHz 

tap input 




Tap Output Amplitude 

FIGURE 3.19 Non Uniformity Data, f _ = 2.5 MHz. 

samp 

R. =10 Kilo-ohms. f. .= 550 KHz 
tap input 
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Tap Output Amplitude 

FIGURE 5 *20 Non Uniformity Data, f =2.5 MHz 

^ samp 

R. - 100 Kilo-ohms. f. .-80 KHz 
tap input 




FIGURE 5.21 Non Uniformity Data, 

®tap“ Kilo-ohms, ^input"* 
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^ap“ ■' 




R. =1 Meg-ohm. f. .=■ 100 KHz 
tap input 
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Input Frequency 30 KHz 100 KHz 
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FIGURE 5.24 Pictorial Presentation of TAD-12 Performance* ^gajup“ ^.5 MHz. ^ Meg-ohm 

V ,* .06 volt peak to peak. Scale; H ® 10 microsec. /div; V « .05 volt/div 
input 



20 KHz 50 KHz 100 KHz I 50 KHz 




A A 



^tap” Kilo-ohms; ^^nput** Volt peak to peak(sinusoidal input) . 
Scale: H*10 raicroaec./div, , V«,2 V/div, 



Tap #1 Output 



Tap #5 Output 




Tap #6 Output 



Tap #11 Output 




FIGURE 3*26 Pictorial Presentation of TAU-12 Performance. 

f “3 MHz; R. =100 Kilo-ohms, 
samp tap 

Vinput* volt peak to peak; ^^nput^* 

Scale: H = 10 microsec. /div, V = .2 volt/div 
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Tap #1 Output 



Tap #4 Output 



Tap #2 Output 



Tap #3 Output 



FIGURE 3.27 





Pictorial Presentation of Reticon TAD-12 



Performance, f = 3*5 MHz. 
samp 

7 .6 Volt peak to peak, 

input 

Scale: H = 10 microsec. /div; 



-100 Kilo-ohms 

tap 



^ input 
V =• .2 



* 50 KHz. 
Volt/div 
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Tap #1 Output 



Tap #2 Output 




Tap #3 Output 



Tap #12 Output 




FIGURE 5»28 Pictorial Presentation of Reticon TAJ)-12, 



Performance 



V, .6 

input 

Scale : H = 



. f - 4 MHz. R, - 
samp tap 

V peak to peak; 
r 10 microaec./div; V = 



100 Kohma. 
50 KHz. 

.2 V/div. 
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Tap #1 Output 



Tap #9 Output 




FIGURE 5.29 



Reticon TAU-12 Performance. 

f = 4.5 MHz. R. = 100 Kilo-ohms, 
samp tap 

Vinput”* ^ peak to peak. ^input°" 

Scale: H = 10 micro sec./div; V * .2 ?/div. 



Tap #1 Output 



Tap ffO. Output 




FIGURE 3 *50 Reticon TAU-12 Performance. 

f =5 MHz. R, =100 Kilo-ohms, 
samp tap 

^input“ ^ ^input’ 50 KHz. 

Scale; H = 10 microsec./div; V = .2 V/div. 
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Input 



Tap #1 Output 




FIGURE 3«31 Reticon TAU-12 Performance. 

f “2.5 MHz. R. “ 10 Kilo-ohms 
samp tap 

V- .6 V peak to peak. f. .=* 

input ^ input 

Upper Scale: H = 5 micro sec. /div; 

Lower Scale: H = 5 microsec. /div; 




FIGURE 3.32 Reticon TAD-12 Performeince 

f = 2.5 MHz. R, =10 Kilo-ohms, 
samp tap 

V. .= .6 V pesUc to T)eak. f. .= 
input ^ ^ input 

Upper Scale: H = .5 microsec. /div; 
Lower Scale: H = .5 microsec. /div; 



100 KHz. 

V = .2 V/div 

V = .1 V/div 



MHz. 

= .2 V/div 
= .1 V/div 
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Input 



Tap #6 Output 




FIGURE 5»33 Reticon TAU-12 Performance. 

f * 2.5 MHz. R. =10 Kilo-ohms, 
samp tap 

Vinput* ^ peak to peak, ^^nput" KHz. 
Upper Scale: H = 1 microsec. /div; V = .2 V/div 
Lower Scale: H = 1 microsec. /div; V = .1 V/div 



Input 



Tap #6 Output 



FIGURE 3»34 Reticon TAI>-12 Performance. 

f = 2.5 MHz; R. =10 Kilo-ohms, 
samp tap 

V. .» .6 V peak to peaik. f, .= 1 MHz 
input ^ ^ input 

Upper Scale: H = .5 microsec. /div; V = .2 ?/div 
Lower Seale: H = .5 microsec. /div; V = .1 v/div 
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Input 



Tap #12 Output 




FIGURE 3.35 



Reticon TAD-12 Performance. 



samp 



2.5 MHz. R, 



tap 



10 Kilo-ohms. 



V. .= .6 V peak to peaUc. f. 
input ^ ^ input 

Upper Scale: H = 5 microsec. /div; V 
Lower Scale: H =• 5 microsec. /div; V 



100 KHz 
» .2 V/div 
- .1 V/div 



Input 



Tap #12 Output 




FIGURE 3.36 Reticon TAD- 12 Performance. 

^samp' ^-5 ■‘tap’ Kilo-ohms. 

’'input" '' 300 KHs 

Upper Scale: H * 2 microsec. /div; V » .2 V/div 
Lower Scale: H = 2 microsec. /div; v = .1 V/div 
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Input 



Tap #3 Output 




FIGURE 3*37 Reticon TAD- 12 Performance. 



f » 3 MHz. R, =10 Kilo-ohms, 
samp tap 

V. ,=> .6 V peak to peak; f. , 
input ^ ^ input 

Upper Scale: H = 2 microsec./div; 
Lower Scale: H = 2 microsec./div; 



= 200 KHz 

V = .5 V/div 

V = .2 V/div 



Input 



Tap #2 Output 




FIGURE 3*30 Reticon TAD- 12 Performance. 

f * 3 * 5 MHz; R, =10 Kilo-ohms, 
samp tap 

V. ,= .6 V peak to peak; f. ,= 200 KHz 
input ^ ^ input 

Scale: H = 2 microsec./div; 7 = .5 V/div 
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Input 



Tap #2 Output 




FIGURE 3*39 Reticon TAD- 12 Performance. 

f » 5 MHz. R. =10 Kilo-ohms 
samp tap 

V, .= .6 V peak to peak. f. .= 200 KHz 
input ^ ^ input 

Upper Scale; H = 5 micro sec./div; 7 = .2 V/div 
Lower Scale; H » 5 microsec. /div; V = .01 V/div 



Input 



Tap #1 Output 



FIGURE 3*40 Reticon TAD-12 Performance. 

f =5 MHz. R. =10 Kilo-ohms 
samp tap 

7. .= .6 7 peak to peak; f. .= 200 KHz 

input input 

Upper Scale; H = 2 microsec. /div; 7 « .2 7/div 
Lower Scale; H = 2 microsec. /div; 7 = .1 7/div 
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Input Frequency 50 KHz 100 KHz 
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FIGURE 5*4^ Second Haxmonic Distortion of Different tap outputs. Signal Input Frequency - 10 KHz. 



20 log(V 2nd harm*/ V fund.) 
- 10 Kilo-ohms 
-100 Kilo-ohms 
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500 KHz £c = 1.25 NHz = 2.5 MHz 
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c. TIME DOMAIN EVALUATION 

The measurement procedure is to use a periodic input 
pulse. Its period is at least twice the total delay. The 
width of the pulse must satisfy the condition: 

Y of clock period £ width of pulse < 1 clock period. 

The pulse amplitude should be within the dynamic range. At 
the output of the summer operational amplifier, one can ob- 
serve twelve output amplitudes. An analog sample and hold 
circuit is used to measure each output amplitude. 

Figures 3-45 through 3-47 illustrate the impulse response 
of TAD- 12 performance at different sampling frequencies: 

fs = 50 KHz, 1.25 MHz, 2.5 MHz 

and tap resistances: 

= 10 Kfi, 100 KQ, 1 Mn 

connected to tap #1, 5, 6 . 

1 . Non- Uniformity 

The non-uniformity data of tap outputs for the impulse 
input is shown in Figures 3-48, 3-49, using 10 Ki^ tap resist- 
ances at different sampling frequencies: 

fs = 25 KHz, 50 KHz, 500 KHz, 1.25 MHz. 

The variation of the tap output amplitudes is the 
following : 



4.21, 


at 


fs - 25 


KHz 


1 . 8 %, 


at 


£3 = 50 


KHz 


4.1%, 


at 


£3 = 500 


KHZ 


10.47% 


, at 


f S = 1 . 


25 MHz. 
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Figures 3-50, 3-51 illustrate the output of 12 taps 
using = 10 sampling frequencies: 

fg = 50 KHz, 1.25 MHz. 

At lower sampling frequency, up to f^ = 1.25 MHz, the 
non-uniformity of tap output is less than 10%; when the samp- 
ling is increased, the uniformity is deteriorated seriously. 
At f^ = 2.5 MHz, the non-uniformity of different taps is 
higher than 20%. 

The interaction between adjacent taps can also affect 
the tap uniformity at high sampling frequencies. Figure 
3-45 (c) , (d) , (ej , (f) shows these effects. For example, 
in Figure 3-45 (dj at f^ = 1.25 MHz, the tap #1 output can 
have an effect on tap #2 before its amplitude reaches the 
zero value. Another example, in Figure 3-45 (f) , at f^ of 
2.5 MHz, the effect of tap #1 output can extend to tap #2 
and tap #3 before its amplitude goes below the fixed pattern 
noise . 

2 . Temporal Noise 

The noise of tap outputs can be grouped into two 
categories : 

- fixed pattern noise 

- random noise. 

The fixed pattern noise is periodic with a period 
related to the sampling frequency used. At f^ of 50 KHz, 
this period is about 6.67 microseconds. At f^ of 500 KHz, 
it is about 80 microseconds, but the number of periods is 
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constant between tap outputs. For example, there are three 
periods of fixed pattern noise between tap #1 output and 
tap #5 output. This is illustrated in Figures 3-45 (a), 

(c) , 3-46 (a), (c) , 3-47 (a), (c) . The peak to peak ampli- 
tude of this noise is about 40 mV for the 10 Kfl load resistors 
and about 25 mV for the 100 Kfl load. 

The random noise varies from one tap to the other. 

It depends also on the clock frequency and tap resistors. 

The random noise affects seriously the tap output, except 
tap #1, at high tap resistor value and high clock frequency. 
For example, using 10 Kfi resistors, the amplitude of random 
noise is 40 mV at clock frequency of f^ = 2.5 MHz and is 
100 mV at f^ = 5 MHz. This is shown in Figure 3-45 (d) , (f) . 
The tap #5 output is deteriorated badly at f^ = 5 MHz for 
R = 100 Kfi and R = 1 Mf2 in Figures 3-46 (d) and 
3-47 (d) . However, the noise on tap #1 is very small at any 
condition. This is illustrated in Figures 3-45 (d) , (f ) , 

3-46 Cd), (f). 

The jitter amplitude observed at the top of tap #1 
output in Figure 3-52 equals to 80 microvolts for a signal 
frequency of 2.83 KHz, f = 100 KHz and for R = 10 Kfi. 

C L 
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Figure 3-5 shows the magnitude frequency response 
of Tap #1, Tap #6 and Tap #8 outputs using 100 Kilo-ohms 
resistors and 1 Meg-ohm resistors. The roll-off starts at 
about 20 KHz input frequency for the = 1 Meg-ohm curve 

and at about 60 KHz input frequency for the R^^p = 100 Kilo- 
ohms curves . 

When the sampling frequency is increased, the fre- 
quency dependence curves are not improved in comparison with 
the theoretical curve. Figures 3-6 and 3-7 illustrate the 
magnitude frequency response of different tap outputs at 
sampling frequency equal to 2.5 MHz using tap resistors: 

R^^p = 10 Kilo-ohms, 100 Kilo-ohms, 1 Meg-ohm. 

With R^^p = 10 Kilo-ohms, the frequency roll-off is 
about 150 KHz. When the signal frequency, f, equals to 200 
KHz, the measured curve is 1 db below the theoretical curve. 
At f = 300 KHz, it is 3 db below the theoretical curve. For 
the R.^ = 1 Meg-ohm curves, the frequency roll-off begins 

about 30 KHz. 

With R^^p = 100 Kilo-ohms, the response curve begins 
to roll off at the input frequency f = 80 KHz. At f = 100 
KHz, the response curve is 2 db below the theoretical curve. 

Figures 3-8 and 3-9 illustrate respectively the 
magnitude frequency response of TAD- 12 output using 10 Kilo- 
ohms tap resistors at sampling frequency f^ equal to 2.75 MHz 
and using 1 Meg-ohm tap resistor at fg equal to 1.25 MHz. 
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50 KHz 1.25 MHz 2.5 MHz 
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Lower Scale: H=1 microsec./div. V=.02 V/div.; (d) Scale: H-1 micros ec./div. V-.1 V/div. 

(e) Upper Scale: H=2 microsec./div. V=1 V/div. , Lower Scale: H=2 microsec./div. V-.02 V/div. 

(f) Upper Scale: H=10 microsec./div. V-.05 V/div,, Lower Scale: H=10 microsec./div. V».1V/div 
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(b) Scale: H»10 mlcrosec./div, V-.5 V/div. ; (c) Upper Scale: H-1 microsec. /div. V«1 V/div. , 

Lower Scale: H=>1 microsec, /div. V»,05 V/div.; (d) Scale: H»5 microsec, /div, V=,02 V/div, | 

(e) Upper Scale: H=2 microsec. /div, V”2 V/div, , Lower Scale: H«2 microsec, /div, V*.01 V/div. ; 

(f) Scale: H»2 microsec. /div. V»,1 V/div, 
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(b) Scale: H*50 micro sec •/div. V».1 V/div. ; (c) Upper Scale: microsec. /div. V»2 V/div. , 

Lower Scale: H-.5 microsec./div. V-.002 V/div. ; (d) Scale; H»5 microsec. /div. V«.01 v/div.; 

(e) Upper Scale; H=.2 microsec. /div. V=1 V/div. , Lower Scale; H=*.2 microsec./div. V“=.002 V/div. » 

(f) Upper Scale; H-10 microsec./div. V-.01 v/div. , Lower Scale; H-10 microsec./div. V=.02 V/div. 








PIGT3HE 3»48 Non Dhifonnity Data. Input Pulse: Amplitude»1 Volt, 
Width=.2 microsec.. Tapping Resistor=10 Kilo-ohms 

(a) P = 25 KHz 

' ^ samp 

(b) P = 50 KHz 

' samp 
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Tap output amplitude 




FIGURE 3»49 Non Uniformity Data* Tapping Resistors“l6 Kilo-ohms 

(a) Input Pulse r Amplitude=1 . 1 Volt, Width=,6 microsec. 



F^^^ -500 KHz 
samp 

(b) Input Pulse; Amplitude-1.1 Volt, 



Width-. 5 raicrosec. 



F 

samp 



1.25 Mz 






8.9 



Output 



Input 



?IGDRE 



Output 



Input 



FIGURE 5 




.50 Impulse Response of Reticon TAD- 12 showing 
12 Tap Outputs. ®-tap* Kilo-ohms. 

Sampling Frequency = 50 KHz. 

Upper Scale: H = 50 microsec. /div; V« .5 V/div 
Lower Scale; H = 50 microsec. /div; 7=1 V/div 




,51 Impulse Response of Reticon TAD-12 Showing 
12 Tap Outputs. Kilo-ohms. 

Sampling Frequency = 1.25 MHz. 

Upper Scale; H = 2 microsec. /div; V = .2 V/div 
Lower Scale; H = 2 micro sec. div/; 7=1 7/div 
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Output 



Tap #1 Jittering 
Waveform 




FIGURE 5»52 Upper Trace: Impulse Response of 12 Tap Outputs 

Scale: H » .1 msec./div; V = ,05 V/div 

Lower Trace: Jittering Waveform of Tap #1 

Scale: H = ,1 msec./div; V » ,02 Y/div 

Sampling Frequency » 50 KHz 
Tapping Resisteuice *10 Kilo-ohms, 
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IV. INVESTIGATION OF FILTER I - PREWHITENING FILTER 



A. THEORY 

Based on the design algorithm described in II.D.l, a 
computer program shown in Appendix A is written to give the 
impulse response of the prewhitening filter for the required 
frequency characteristics. Using this program, a prewhitening 
filter is designed for the following conditions: 

f^ = 2.1 KHz f = 1/2 f = 625 KHz 

o p s 

R = 0.1 Slope A = 1 

Length of filter: M = 11 taps 




The tapping coefficients improved by Hamming window are 
listed in Table IV. I. 

The magnitude frequency response of the theoretical pre- 
whitening filter is shown in Figure 4-1. It is plotted in 
linear scale on Y-axis versus normalized frequencies, 
f/fg on X-axis. 
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TABLE IV. I Theoretical Tapping Coefficients § Resistors 



k 

-5 

-4 

-3 

-2 

-1 

0 

1 

2 

3 

4 

5 



-0.00091 

0.0 

-0.01219 

0.0 

-0.21407 

1 . 

-0.21427 

0.0 

-0.01219 

0.0 

-0.00091 



(Kilo-ohms) 

-16478.516 

open 

-1225.517 

open 

-65.005 

10 . 

-65.005 

open 

-1225.517 

open 

-16478.516 



A, Measured 
k 

0.02564 

0.15385 

0. 05128 

0 . 0 

-0.23076 

1 . 

-0.19231 
0.02564 
0.02564 
0. 02566 
-0.07692 



B. EXPERIMENT 



1 . Tapping Resistors According to Theoretical 
Design with ^min - 10 Kilo-ohms 

A minimum value of 10 Kilo-ohms is set for the tapping 



resistances. Appendix F indicates the following expression 
for Rj^, 



R, = (lOK + r) 0 



^max 



) - r 



where r is the internal output resistance. For the TAD-12, 
r is about 5 Kilo-ohms. The tapping resistors corresponding 
to the designed prewhitening filter are given in Table IV. 1 
and used to implement the filter using the 11 taps of 
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f 




Reticon TAD-12. Since their theoretical coefficients 
for k higher than 5 have negligible values, the implementa- 
tion of the filter on a longer delay line, number of taps 
higher than 11, is not necessary. 

The measured impulse response of this filter is 
illustrated in Figure 4-2 at three different sampling fre- 
quencies : 

fg = 50 KHz, 500 KHz, 1.25 MHz 

The pulse inputs used have the amplitude and width listed 
in Table IV. II. 

TABLE IV. II Pulse Inputs 
fg Amplitude Width 

50 KHz 1.5V 8 micros ec. 

500 KHz 1.5 V .8 microsec. 

1.25 MHz 1.5 V .3 microsec. 

Using the sample and hold circuit, the impulse 
response, A^^ is measured and the result is listed in Table 
IV. 1. The Aj^ measured and the theoretical Aj^ are quite 
different for some taps. The maximum deviation occurs be- 
tween tap #7 and #8 outputs. Its peak amplitude is about 
301 of the output at tap #6. 

The theoretical magnitude frequency response and 
measured response curves for three different sampling fre- 
quencies are plotted on Figure 4-3. The response curves for 
f^ of 1.25 MHz and 500 KHz show some wiggling shape due 
probably to the non-uniformity of tap output and the clock 
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noise which cause the unsymmetrical tapping coefficients. 

The response curve for f equal to 50 KHz is rather satis- 
factory if it is multiplied by a factor of 1.07. 

2 . Tapping Resistors Adjusted to Yield 

Better Impulse Response. Rmin = 10 Kilo-ohms 

Since only three taps have strong contribution to 

the filter, the adjustment is done only on the taps #5, #6, 

#7 by using the sample and hold circuit. The other taps are 

left open. The tapping resistors used are listed in Table 

IV. III. 



TABLE IV. I II. Tapping Resistors 

k A. (Kilo-ohms) 

^ fs = 50 KHz fs = 500 KHz fs = 1.25 MHz 



1 


-0.21427 


28 


50 


52 


0 


1. 


10 


10 


10 


1 


-0.21427 


28 


35 


100 



The amplitude and width of inputs pulse are listed in 
Table IV. II. 

The impulse response of this filter is illustrated in 
Figure 4-4, at three different sampling frequencies. 

The theoretical and measured curves of prewhitening 
filter magnitude frequency response are shown in Figure 4-5. 
The response curves for f^ equal to 50 KHz, 500 KHz and 1.25 MHz 
still show some wiggling. The gain of these curves are lower 
than the last curves in Figure 4-3. If these three curves 
are multiplied by a factor of 1.22, they are somewhat satis- 
factory in comparison with the theoretical curve. 
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3 . Tapping Resistors Adjusted to Yield 

Better Impulse Response. Rmin = 50 ohms 

A minimum value of 50 ohms is set for the tap #6. 

The taps #5 and #7 are adjusted by using sample and hold 

circuit, to yield better impulse response. The other taps 

are left open. The tapping resistors used are listed in 

Table IV. IV. 

TABLE IV. IV. Tapping Resistors 

k A]j (Kilo-ohms) 

fg = 50 KHz fg = 500 KHz fg = 1.25 MHz 



1 


-0.21427 


00 


2.8 


00 


0 


1. 


0.05 


0.05 


0.05 


1 


-0.21427 


2.8 


1.8 


11.5 



The pulse inputs used have the amplitudes and widths 
listed in Table IV. II. 

Figure 4-6 illustrates the impulse response of the 
filter at sampling frequencies; 

fg = 50 KHz, 500 KHz, 1.25 MHz 

The theoretical and measured curves of prewhitening 
filter magnitude frequency response are shown in Figure 4-7. 
The wiggling of the three measured curves is much reduced in 
comparison with the previous curves. If the experimental 
curve f^=l*25 MHz is multiplied by a factor of 1.2 and curves 
fg = 50 KHz, 500 KHz multiplied by a factor of 1.35, the 
agreement is somewhat satisfactory. 
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4 . Tapping Resistors Obtained by Experimental 

Cut § Try Based on Frequency Response Adjustment 

The procedure consists of the use of a chirped wave- 
form from DC to 1 MHz as the input of the filter, and adjust 
the displayed frequency response to get the minimum wiggling 
by adjusting the tapping resistances. The result obtained 
is shown in Figure 4-9. The wiggling of the three response 
curves is much reduced. Raising these curves by a factor of 
1.27 can improve the agreement between the theoretical and 
measured curves. 

The impulse response corresponding to these curves 
is unsymmetrical and has the values listed in Table IV. V. 

TABLE IV. V Tapping Coefficients 
k Aj. Rk (Kilo-ohms) 

-1 -0.1694 24.812 

0 1. .05 

1 -1 .05 

The amplitude and width of pulse inputs used are 
listed in Table IV. II. 
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COL dCZ OO'i 006 OUi CIO 0X0 014 




FIGURE 4.1 Theoretical Amplitude Frequency Response of Prewhitening 
Filter. Length =* 11. Hamming window. 

X - Scale “ .1 units/inch for normalized frequency 

Y. - Scale » .2 units/inch for Amplitude frequency response. 
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(c) Upper Scale: H-.3 miorosec./div. V»1,5 V/div.; Lower Scale: H»5 microsec. /div. V»2 v/di 
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FIGURE 4.3 Prewhitening Filter Magnitude Frequency Response. 

Tapping Resistors According to Theoretical Design, Kilo-ohms 



50 KHz 500 KHz 1.25 MHz 
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(c) Upper Scale: H»5 microsec./div. V*,02 V/div.; Lower Scale: H*5 mlcroseo./div, V«1 V/div 
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FIGURE 4/. 5 Prewhitening Filter Magnitude Frequency Response. 

Tapping Resistors Adjusted to Yield Better Impulse Response. R^^^ - 10 Kilo-ohms 
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(c) Upper Scale: H=*5 microaec*/div* V«»5 V/div.; Lower Scale: H»5 laioroaec./div. V®2 V/div 
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FIGURE ^-7 Prewhitening Filter Magnitude Frequency Response. 
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(c) Scalet H*»5 mlcrosec./div. V*»1 V/div, 




H<f) 
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FIGURE i$.9 Prewhitening Filter Magnitude Frequency Response. 

Tapping Resistors Obtained by Cut and Try Based on Frequency Domain Response. 



V. INVESTIGATION OF FILTER II - DEWHITENING FILTER 



A. THEORY 

The design algorithm of the dewhitening filter is 
described in II. D. 2. The computer program in Appendix B is 
used to calculate the impulse response of the filter from 
the given frequency response. A dewhitening filter is 
designed, based on the impulse response of the prewhitening 
filter of paragraph IV. A. The frequency characteristics 
are the following: 

fo = 2.1 KHz " I = 625 KHz 

Length of filter M = 11 taps 




frt f f normalized 

0 p 

The tapping coefficients are improved by a Hamming window 
and by a factor of 1.13 and are listed in Table V.I. 

The magnitude frequency response of the theoretical pre- 
whitening filter is shown in Figure 5.1. It is plotted in 
linear scale on Y-axis versus normalized frequencies, 
f/fs on X-axis. 
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TABLE V . I . Theoretical Tapping Coefficients § Resistors 



k theor. (Kilo-ohms) (Kilo-ohms) 

with Rjnin “ K ohms, with Rmin = 50 ohms 



-5 


0.00113 


13269.336 


4464.026 


-4 


0.00632 


2368.417 


794.050 


-3 


0.01156 


1292.577 


431.851 


-2 


0.04846 


304.661 


99.252 


-1 


0.19681 


71.215 


20,659 


0 


1. 


10. 


. 05 


1 


0.19681 


71.215 


20.659 


2 


0.04844 


304.661 


99.252 


3 


0.01156 


1292.577 


431.851 


4 


0.00632 


2368.417 


794.050 


5 


0.00113 


13269.336 


4464.026 


EXPERIMENT 






1. 


Tapping Resistors 


According to 


Theoretical 




Design with Rmin 


= 10 Kilo-ohms 





A minimum value of 10 Kilo-ohms is set for the 
tapping resistors. Using the expression of Rj^ given in 
Appendix F, the tapping resistances are calculated and 
listed in Table V.I in the third column. They are used to 
implement the dewhitening filter using 11 taps of Reticon 
TAD- 12. 

The impulse response of this filter is illustrated 
in Figure 5-2 at three different sampling frequencies: 

fs = 50 KHz, 500 KHz, 1.25 MHz. 
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TABLE V.II Theoretical ^ Measured Coefficients 

\ 

k A]^ theor. A]^ Measured 

fs = 50 KHz fs = 500 KHz fg = 1. 25 MHz 



-5 


0.00113 


0.16216 


0.10869 


0.11667 


-4 


0.00632 


0.16216 


0.17391 


0.20830 


-3 


0.01156 


0.16216 


0.15217 


0.25 


-2 


0.04844 


0.2027 


0.19565 


0.2708 


-1 


0.19681 


0.39119 


0.41304 


0.41667 


0 


1. 


1. 


1. 


1. 


1 


0.19681 


0.39189 


0.39130 


0.41667 


2 


0.04844 


0.18918 


0.13043 


0.0 


3 


0.01156 


0.16216 


0.13043 


0.16667 


4 


0.00632 


0.17567 


0.13043 


0.16667 


5 


0.00113 


0.16216 


0.10869 


0.16667 


The impulse 


inputs are noted in Table 


IV. II for 


different 


sampling 


frequencies . 


The tapping 


coefficients 


are measured with the sample and 


hold circuit 


and the result 


are listed 


. in Table 


V.II. The theoretical and measured Aj^ 


are quite 


different 


for some taps 


due to the 


fixed pattern 


noise and 


non-uniformity of tap outputs. The 


maximum noise 


amplitude 


is observed before the 


tap #1 output and has 



amplitude of about 30% of Aq. 

The theoretical magnitude frequency response and 
measured response for three different sampling frequencies 
are plotted in Figure 5-3. The response curves for f^ 
of 500 KHz and 1.25 MHz are somewhat in agreement with the 
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theoretical curve. Their difference is about 12% for f 
less than 0.2 f^ and it becomes worse than 40% for input 
signal frequency higher than 0.25 f^. For f^ of 50 KHz, 
the result is quite satisfactory. 

2 . Tapping Resistors Adjusted to Yield Better 

Impulse Response. Rmin = 10 Kilo-ohms 

Since the coefficients, for k greater than 3, 
have negligible values in comparison with the noise, the 
adjustment is done only on seven tap outputs. The adjust- 
ment is quite delicate for the tap outputs which are greater 
than or equal to 5, because of the interaction between them. 
Experimental tap resistances are listed in Table V.III. A 
slight modification has to be done on the values of some 
tap resistors in order to get better impulse responses for 
different sampling frequencies. The sample and hold circuit 
is used to make the adjustment. 

The amplitude and width of inputs pulse used are 
listed in Table IV. II. 

The impulse response of this filter is illustrated 
in Figure 5-4 at three different sampling frequencies. The 
maximum noise occurs near tap #1 output; its peak amplitude 
is about 19% of A^. 

The theoretical and measured curves of dewhitening 

filter magnitude frequency response are shown in Figure 5-5. 

The wiggling of the response curves is reduced in comparison 

with the curves in Figure 5-3, but the difference between the 

theoretical and measured curves (f =1.25 MHz and 500 KHz) 

'■ s 
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is still high, greater than 30%, for the input signal 
frequency f higher than 0.25 f ^ . At low sampling frequency, 
50 KHz, the result is quite satisfactory. 

TABLE V.III. Adjusted Coefficients § Resistors 



k Ak Rk 

f =50 KHz f =500 KHz f =1.25MHz f =50 KHz f^ = 500 KHz f =1 . 25MHz 
s s s s s s 



-5 


0.09836 


0.02083 


0.125 


13269. 


13269. 


13269. 


-4 


0.10655 


0.02083 


0.1 


2368. 


2368. 


2368. 


-3 


0.10655 


0.04167 


0.075 


1292.5 


1292.5 


1292.5 


-2 


0.09836 


0.0625 


0.075 


300. 


300. 


300. 


-1 


0.18852 


0.20833 


0.175 


110. 


109. 


110. 


0 


1. 


1. 


1. 


5. 


9.5 


5. 


1 


0.18852 


0.20833 


0.2 


110 


110. 


110. 


2 


0.09836 


0.0625 


0.1 


311. 


311. 


320. 


3 


0.10655 


0.02083 


0.075 


1292.5 


1295.5 


1292.5 


4 


0.11675 


0.02083 


0.075 


2368. 


2368. 


2368. 


5 


0.10655 


0.02083 


0.075 


13269. 


13269. 


13269. 




3. Tapping 


Resistors According 


to Theoretical Design. 






Rmin 


50 ohms. 












A minimum value 


of 50 ohms 


is set for 


the tap #6. 


Using 


the 


formula described in 


Appendix F 


, one can 


calculate the 



theoretical tapping resistances whose values are listed in the 
fourth column of Table V.I. 

The input pulses have the amplitude and width listed 
in Table IV. II. Figure 5-6 illustrates the impulse response 
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of the filter at sampling frequencies: 



fg = 50 KHz, 500 KHz, 1.25 MHz 
The measured coefficients are listed in Table V.IV: 



TABLE 


V.IV. Measured 


Coefficients (R„,- 
^ min 


=50 ohms) 


k 


f =50 KHz 
s 


^ -^k 

fg=500 KHz 


f^=1.25 MHz 
s 


-5 


0.11111 


0.06667 


0.22826 


-4 


0.11111 


0.06667 


0.08695 


-3 


0.11111 


0.06667 


0.09782 


-2 


0.16667 


0.13333 


0.15217 


-1 


0.27778 


0.26667 


0.34782 


0 


1. 


1. 


1. 


1 


0.22223 


0.26667 


0.48913 


2 


0.13888 


0.13333 


0.20652 


3 


0.11111 


0.06667 


0.13043 


4 


0.11111 


0.06667 


0.04347 


5 


0.11111 


0.06667 


0.07608 



The maximum peak amplitude of noise measured is about 
131 of for f^ equal to 1.25 MHz and 500 KHz and 19% of A^ 
for f^ equal to 100 KHz. 

The theoretical and measured curves of dewhitening 
filter magnitude frequency response are shown in Figure 5-7. 
The response curve for f^ equal to 50 KHz is similar to the 
theoretical one. The curves f^ equal to 1.25 MHz and 500 KHz 
have some wiggling effect for input frequencies higher than 



.29 f3. 
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000 000 OiO OiS 020 02b 030 




FIGURE 5.1 Amplitude Frequency Response of Dewhitening Filter. 
Length * 11 / Hamming window 

X - Scale = .1 units/inch for normalized frequency 
Y - Scale * .5 units/inch for amplitude frequency 
response. 
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(c) Upper Scale; H»5 micro sec. /div. V-,01 V/div. ; Lower Scale: H-5 microsec. /div, V«1 V/di 
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FIGURE 5.3 Dewhitening Filter Frequency Response. 

Tapping Resistors According to Design. “10 Kilo-ohms 



50 KHz 5D0 KHz 1,25 MHz 
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(c) Upper Scale: H-5 microaec./div, V-.01 V/div, ; Lower Scale: H»5 micro sec. /div. V»2 V/di 
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FIGURE 5.5 Dewhitening Filter Magnitude Frequency Response. 

Tapping Resistors Adjusted to Yield Better Impulse Response . - 10 Kilo-ohms 
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(c) Upper Scale: H=*5 microsec. /dlv, V».02 V/div.; Lower Scale: H-5 microsec. /div. V-1 V/di 
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VI. INVESTIGATION OF FILTER III - BANDPASS FILTER 



A. THEORY 

Based upon the design algorithm described in II. D. 3, a 
computer program shown in Appendix C is written to calculate 
the impulse response of the bandpass filter based on the 
required frequency response. Using this program, a bandpass 
filter is designed for the following conditions: 

= .1996 £2 = .2004 

Center frequency f ^ = . 2 Q = 250 

Amplitude A = 1. 

Number of taps M = 11. 




The tapping coefficients improved by Hamming window are 
listed in Table VI. I. 

The magnitude frequency response of the theoretical band- 
pass filter is shown in Figure 6-1, in linear scale, and in 
Figure 6-2, in logarithmic scale versus normalized frequencies, 

f/f. . 
s 
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TABLE VI. I. Theoretical Coefficients and Resistors 



k 




= k) Rv(R 
k min k^ min 


-5 


0.09863 


147.083 


46.201 


-4 


0.07378 


198.307 


63.447 


-3 


-0.38389 


34.074 


8.155 


-2 


-0.59148 


20.360 


3.538 


-1 


0.28647 


47.362 


12.628 


0 


1. 


10. 


.05 


1 


0.28647 


47.362 


12.628 


2 


-0.59148 


20.360 


3. 538 


3 


-0.38389 


34.074 


8.155 


4 


0.07378 


198.307 


63.447 


5 


0.09863 


147.083 


46.201 


EXPERIMENT 






1. 


Tapping Resistors 


According to Theoretical 




Design with Rmin ' 


= 10 K ohms. 





A minimum value of 10 Kilo-ohms is set for the 
tapping resistors. The expression of given in Appendix F 
is used to calculate the tapping resistances. These values 
are listed in Table VI . I in the third column and are used to 
implement the bandpass filter using the 11 taps of Reticon 
TAD- 12. 

The impulse response of this filter is illustrated 
in Figure 6-3 at three different sampling frequencies: 

fg = 50 KHz, 500 KHz, 1.25 MHz 
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TABLE VI. II. Theoretical § Measured Coefficients 



k 


A-^ Theor. 


f =50 KHz 
s 


A]q Measured 

f =500 KHz 
s 


f =1.25 MH 
s 


-5 


0.09863 


0.34444 


0.30508 


0.34483 


-4 


0.07378 


0.33333 


0.209339 


0.06896 


-3 


-0.38389 


-0.5 


-0.66102 


-0.86306 


-2 


-0.59148 


-0.6 


-0.66102 


-0.51724 


-1 


0.28647 


0.48888 


0.49152 


0. 72413 


0 


1. 


1. 


1. 


0.86206 


1 


0.28647 


0.44444 


0.54237 


0.4 


2 


-0.59148 


-0.65555 


-0.76271 


-1. 


3 


-0.38389 


-0.48888 


-0.71186 


-0.34482 


4 


0.07378 


0.31111 


0.13559 


0.51724 


5 


0.09863 


0 . 26666 


0.13559 


0.24138 




The input 


pulses used 


are listed in Table 


IV. II for 


different sampling 


frequencies 


. The tapping coefficients 



measured with the sample and hold circuit are shown in Table 
VI. II. The theoretical and measured are quite different 
for some taps due to the fixed pattern noise and the non- 
uniformity of the tap outputs. The maximum noise observed 
occurs between tap #3 and #4 outputs. Its peak amplitude is 
about 30% of A^. 

The theoretical magnitude frequency response and 
measured response curves for three different sampling fre- 
quencies are plotted in Figure 6-4. The main lobes of measured 
response curves are rather similar to the theoretical one. 
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But there are some differences at the side lobes which have 
their peak magnitudes as follows: 

- 15% of maximum amplitude for f^ = 50 KHz 

- 35% of maximum amplitude for f^ = 500 KHz 

- 51% of maximum amplitude for f^ = 1.25 MHz. 

There are also some wiggles for input frequencies less than 

.15 f^. 
s 

2 . Tapping Resistors Adjusted to Yield Better 

Impulse Response, Rmin = 10 Kilo-ohms 

Since for k greater than 3, the coefficients have 
negligible values in comparison with the noise, the adjustment 
is done only on 7 tap outputs. Obtaining the exact values of 
the theoretical design is very delicate because of the inter- 
action between tap outputs. The adjusted coefficients and 
the corresponding tap resistances are listed in Table VI. III. 
Some modification has to be done on the values of some tap 
resistors in order to get better impulse response for differ- 
ent sampling frequencies. 

The amplitude and width of input pulses are listed 
in Table IV. II. 

The impulse response of this filter is illustrated 
in Figure 6-5 at three different typical sampling frequencies. 
The maximum noise observed has its peak amplitude about 19% 
of Aq. 

The theoretical and measured curves of bandpass filter 
magnitude frequency response are shown in Figure 6-6. The 
wiggling of the frequency response at input frequencies less 
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than .15 is much reduced for = 50 KHz, but the side 
lobe is still high at 15% of the maximum amplitude. For = 
1.25 MHz and £^ = 500 KHz, the peak amplitudes of side lobes 
are increased to 29% of maximum amplitude. 

TABLE VI. III. Adjusted Coefficients ^ Resistors CRmin~ 10 



k A^ Rk (Kilo-ohms) 

£^=50KHz £^=500KHz £^=1.25MHz £ =50KHZ £ =500KHz £^=1.25MHz 

s s s s s s 



-5 


.07246 


.01587 


.07317 


147. 


147. 


147. 


-4 


.08695 


.04762 


.04878 


198. 


198. 


198. 


-3 


-.20289 


- . 55555 


-.34146 


48.07 


40. 


60.11 


-2 


-.62318 


- .62492 


- . 53658 


23. 


20. 


16. 


-1 


.28985 


.30156 


.26829 


47. 


75. 


210. 


0 


1. 


1. 


1. 


8. 


8. 


6. 


1 


.28985 


.30156 


.26829 


47. 


75. 


9.8 


2 


- . 62318 


- .63492 


- . 53658 


23. 


24.8 


8.5 


3 


-.20289 


- . 55555 


-.34146 


48.07 


40. 


23. 


4 


. 08695 


.04762 


. 07317 


198. 


198. 


198. 


5 


.07246 


.01587 


- .04878 


147. 


147. 


147. 



3 . Tapping Resistors According to Theoretical 

Design. Rmin = 50 ohms 

A minimum value of 50 ohms is set for the tapping 
resistance. Using the formula described in Appendix F, the 
tapping resistances are calculated and have the values listed 
in Table VI . I , fourth column. 

The pulse inputs used have the amplitude and width 
listed in Table IV. II. Figure 6-7 illustrates the impulse 
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the impulse response of the filter at sampling frequencies: 



f = 50 KHz, 500 KHz, 1.25 MHZ 
s 



measured 


coefficients Aj^ 


are listed 


in Table VI. IV. 


TABLE VI 


.IV. Measured 


Coefficients 


(R 50 ohms) 

min 


k 


fg=50 KHz 


■^k 

f =500 KHz 
s 


f =1.25 MHz 
s 


-5 


.03211 


.03137 


. 82810 


-4 


.05504 


. 07058 


.82062 


-3 


- .37155 


- . 27451 


. 76083 


-2 


- . 58256 


-.40392 


. 74439 


-1 


. 24312 


.22745 


.84603 


0 


1. 


1. 


1. 


1 


. 24312 


.24314 


. 84603 


2 


-.58256 


-.34509 


.73243 


3 


- .37155 


-.38823 


. 73094 


4 


.05504 


.04313 


.83408 


5 


.05504 


.09803 


.81464 



The maximum peak amplitude of noise measured is about 18% of Aq. 

The theoretical and measured curves of bandpass filter 
magnitude frequency response are shown in Figure 6-8. The 
wiggling effect of the response curves are strongly diminished 
in comparison with previous curves in Figures 6-4 8 6-6. The 
amplitudes of side lobes are also reduced and have the follow- 
ing values : 
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50 KHz 



7.5% of maximum amplitude for f^ = 

- 11% of maximum amplitudes for f^ = 500 KHz 

- 12.5% of maximum amplitude for f^ = 1.25 MHz. 

4 . Tapping Resistor Adjusted to Yield Better 
Impulse Response. Rmin = 50 ohms 

a. Driving Circuit Using RCA/CMOS Flip-Flop 

A minimum value of 50 ohms is set for tap #6. 

The tap outputs are adjusted by sample and hold circuit to 

yield better impulse response. The tapping coefficients and 

resistances adjusted are listed in Table VI. V. 

TABLE VI. V. Adjusted Tapping Coefficients § Resistors 



k Ak Rv (Kilo-ohms) 

f =50KHz f =500KHz f =1.25 MHz f =50KHz f =500KHz f =1.25 MHz 
s s s s s s 



-5 


.16733 


. 08532 


.01234 


46.2 


46.2 


46.2 


-4 


. 06374 


.07508 


. 03086 


63.45 


63.45 


63.45 


-3 


- .37450 


- .33673 


-.25308 


6.25 


4.05 


4.05 


-2 


- .60557 


-.59183 


-.59876 


1.2 


1.2 


. 5 


-1 


. 28287 


.23849 


.22222 


23.023 


25.023 


10. 


0 


1. 


1. 


1. 


. 05 


.05 


. 05 


1 


. 28287 


.23849 


.23457 


6.1 


6.1 


4.1 


2 


-.60557 


-.59183 


-.59876 


1.083 


1.083 


1.08 


3 


-.37450 


-.33673 


- .23457 


4. 083 


4.55 


1.083 


4 


.06374 


.07508 


. 04321 


63.45 


63.45 


63.45 


5 


.16733 


.08532 


.02469 


46.2 


46.2 


46.2 


The 


input pulse used 


have the 


amplitudes and 


widths 


listed in 



Table IV. II. Figure 6-9 illustrates the impulse response of 
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the filter at three different typical sampling frequencies. 
The maximum noise has its peak amplitude about 15% of A^. 

The theoretical and measured curves of bandpass 
filter magnitude frequency response are shown in Figure 6-10. 
The wiggling does not appear at the main lobe but the ampli- 
tudes of the side lobes are s.till high: 



- 15% 


of maximum 


amplitude 


for 


£3 = 500 


KHz 


- 10 % 


of maximum 


amplitude 


for 


£3 = 500 


KHz § 1.25 MHz 



b. Driving Circuit Using INSELEX CMOS/SOS Flip-Flop 

A new circuit with faster flip-flop is made and a 

minimum value of 50 ohms is set for tapping resistance. The 

adjustment is done for f^ equal to 500 KHz; the coefficients 

and resistances obtained are used also for f of 50 KHz and 

s 

1.25 MHz. Table VI. VI lists these coefficients and resistances. 

TABLE VI. VI Adjusted Coefficients ^ Resistances 
k A]^ R]^ (Kilo- ohms) 



-5 


. 03536 


46.2 


-4 


.01571 


63.45 


-3 


- .40078 


3.8 


-2 


-.61689 


1.56 


-1 


.29194 


9.023 


0 


1 . 


.05 


1 


.29273 


2.32 


2 


- . 61297 


1.263 


3 


- .40078 


5.000 


4 


.06287 


63.45 


5 


.07269 


46. 2 
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The pulse inputs listed in Table IV. II are used. Figure 6-11 
illustrates the impulse response of the filter at sampling 
frequencies . 

f = 50 KHz, 500 KHz, 1.25 MHz. 
s * * 

The maximum peak amplitude noise observed is about 21.6% of 
^o* 

The theoretical and measured curves of bandpass 
filter magnitude frequency response are shown in Figure 6-12. 
There is no wiggling at the main lobe region and some improve- 
ment is obtained at the side lobes, which have the following 
amplitudes : 

6% of maximum amplitude for f^ = 1.25 MHz 

- 10% of maximum amplitude for f^ = 500 KHz 

- 14% of maximum amplitude for f^ = 50 KHz. 

Using the computer program shown in Appendix E, 

the measured magnitude frequency response for f^ equal to 
500 KHz and the response corresponding to the measured impulse 
response listed in Table VI. VI are plotted in Figure 6-13 with 
DB scale on Y axis and normalized frequencies on X axis. The 
main lobes of the two curves are the same but there are some 
differences at the side lobe regions. This shows that the 
technique using the sample and hold circuit to adjust the 
tapping resistor to get better impulse response is quite 
accurate and very useful to the filter implementation. 
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5 . Distortion Measurements of Bandpass Filter 

For the case described in V.B.3 where the tapping 

resistances are set according to theoretical design with 

R . equal to 50 ohms, the second harmonic distortion de- 
min ^ ’ 

fined by: 



DB 



20 1 P 2nd harmonic voltage 
fundamental voltage 



are illustrated respectively in Figures 6-14 and 6-15, for 
sampling frequency of 500 KHz and three different peak to 
peak input voltages: 



V. = 3V, 0.8V, O.IV 
in * ’ 

From the Figure 6-14 where the input frequency is equal to 
50 KHz, the second harmonic and clock noise at 85 KHz have 
the relative values expressed in percentage of fundamental 
as shown in Table VI. VI I. 



TABLE VI. VI I Second Harmonic 8 Clock Noise 



V. ppk 
in 


3 Volts 


. 8 Volts 


.1 Volt 


2nd Harmonic 


16.7% 


6.5% 


6. 8% 


Clock Noise 


5.8% 


14.9% 


88.8% 



When the amplitude of the input voltage decreases, the rela- 



tive value of second harmonic decreases but the relative value 



of clock noise increases. In the Figure 6-14c, one can observe 



many spurious responses which occur at 63 KHz, 85 KHz, 112 KHz, 



126 KHz, and are due principally to the clock noise and fixed 



pattern noise. 
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oeo S2C 020 OTC CTO 5C0 




FIGURE 6.1 Amplitude Frequency Response of Bandpass Filter 
Length = 11 / Hamming Window 

X - Scale = .2 units/inch for normalized frequency 
Y - Scale = 10 units/inch for amplitude frequency 
response. 
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FIGURE 6*2 Amplitude Frequency Response of Bandpass Filter . 

Len^h = 11 / Hamming Window 

X - Scale » .2 units/inch for normalized frequency 
Y — Scale » 10 units/inch in DB for amplitude frequency 
response. 
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FIGURE ^,4 Beindpass Filter Magnitude Frequency Response, 

Tapping Resistors According to Design, R . =10 Kilo-ohms 



50 KHz 500 KHz 1.25 MHz 
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FIGURE 6,6 Bandpass Filter Magnitude Frequency Response, 

Tapping Resistors Adjusted to Yield Better Impulse Response, R ^^ ^ - 10 Kilo-ohms 



50 KHz 500 KHz 1.25 1fflz 
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FIGURE 6.8 Bandpass Filter Magnitude Frequency Response. 

Tapping Resistors According to Design. R “ 50 ohms 



50 KHz -500 KHz 1.25 MHz 
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FIGURE 6.10 Bandpass Filter Habitude Frequency Response. 

Tapping Resistors Adjusted to Yield Better Impulse Response. R -50 ohms 

min 
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.FIGURE 6*12 ijandpass Filter Magnitude Frequency Response. 

Tapping Resistors adjusted to Held Oetter Impulse Response. R^^ = 50 



-oao -03^ -030 -02& -020 -015 -OlO -005 




FIGURE 6.13 



Bandpass Filter Magnitude Frequency Response. 

Tapping Resistors Adjusted Using Impulse Response. 

Rmin = 50 ohms. Q = ,250,. J 
X^*^-^ScaIe = .1 units/inch £c 

Y - Scale = 5 units/inch in DB 
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FIGURE 6.15 Second Harmonic Distortion of Bandpass Filter with Tapping Resistors According 

to Design and 50 ohms for Different Input Signal Amplitudes. fggj„p - 500 KHz 



6. Observations 



The characteristic feature of the filter realized 
by the electronic tapped analog delay line is that one can 
scan the center frequency of the bandpass filter to the left 
or right by simply decreasing or increasing the clock fre- 
quency. Figure 6-16 illustrates this feature for Q = constant 
at two different sampling frequencies: 




FIGURE 6-16. Bandpass Filters Amplitude Frequency 
response at 2 different sampling 
frequencies f =500 KHz § f =835 KHz. 
Scale : H = . 2^nsec/Div V ^ . 1 V/Div. 



Another observation is that the transversal filter 
implementation is not very sensitive to high specifications 
for moderate number of taps. The tapping coefficients are 
almost the same for different values of Q required. The 
impulse response of bandpass filters for: 

f^ = 0.2 fg (f^ = center frequency) 

Q = 50, 250, 500 
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using 11 taps is listed in Table VI. VIII, and using 23 taps 
in Table I. IX. 



TABLE VI. VI I I Tapping Coefficients (11 Taps) 

k Ak 

Q=50 Q=250 Q=500 



0 


1. 


1. 


1. 


1 


0.28644 


0.28647 


0.28640 


2 


-0.59140 


-0.59148 


-0.59148 


3 


-0.38382 


-0.38389 


-0.38395 


4 


0.07375 


0.7378 


0.07378 


5 


0.09857 


0.09863 


0.09863 


TABLE 


VI. IX Tapping 


Coefficients 


(23 Taps) 


k 


Q=50 


Ak 

Q=250 


Q=500 


0 


1. 


1. 


1. 


1 


0.30373 


0.30377* 


0.30370 


2 


-0.75475 


-0.75486 


-0.75486 


3 


-0.69072 


-0.69086 


-0.69095 

\ 


4 


0.23217 


0.23227 


0.23227 


5 


0.63317 


0.63357 


0.63358 


6 


0.15702 


0.15717 


0.15718 


7 


-0.31184 


-0.31222 


-0.31224 


8 


-0.22189 


-0.22225 


-0.22226 


9 


0.05648 


0.05660 


0.05660 


10 


0.11777 


0.11807 


0.11808 


11 


0.02596 


0.02604 


0.02605 
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VI. CONCLUSION 



Digital non- recursive filters have been extensively 
developed for the implementation of "finite impulse response" 
(FIR) filters using software on digital computers. In recent 
years, new CTD (charge transfer device) tapped delay lines 
make it possible to implement real time FIR filters using these 
hardware devices. However, because the signal is sampled 
analog instead of digital and because the device principle 
and technology are new, careful investigation is needed for 
the successful development of these hardware sampled analog 
signal processors. In this thesis, the new Reticon TAD-12 
tapped delay line with external programmable tap weights is 
investigated in detail. In addition to a thorough evaluation 
of its device performance, the properties of three sampled 
analog FIR filters using TAD- 12 are investigated. They are 
the prewhitening, dewhitening and bandpass filters. 

Both the frequency response of TAD- 12 devices have been 
evaluated for sampling frequency up to 2.5 MHz. It was found 
that the output showed the following limitations: 

1. Frequency response limitation: 

Tap outputs generally roll-off faster than predicted 
by the effect of sample and hold. The roll-off depends on 
the location of the tap and its tapping resistance. 

2. Time response limitation: 
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Output at one tap generally takes more than the delay 
period between adjacent taps to settle down and affect the 
output of the following taps. The settling time depends on 
the location of the tap, its tapping resistance and the 
sampling frequency. 

3. Non-uniformities: 

Tap outputs are not uniform among the taps. The non- 
uniformity varies with the location of the tap, its tapping 
resistance and signal frequency. 

4. Loading effect: 

Tap outputs vary with tapping resistance. However, 
it should be pointed out that these limitations are large at 
higher sampling frequencies and signal frequencies. At 
sampling frequencies less than 100 KHz and signal frequencies 
less than 10 KHz, TAD-12 behaves satisfactorily. 

The properties of three sampled analog FIR filters using 
TAD-12 have been evaluated following three design procedures. 

1. Tap resistances are selected using the same design 
procedure developed for digital non-recursive filters. 

2. Tap resistances are corrected using the frequency response. 

3. Tap resistances are corrected using the impulse time 
response . 

It was found that the device limitation causes considerable 
deviation of the filter performance from theoretical prediction 
especially when the sampling frequency is high. Furthermore, 
the deviation is relatively smaller in bandpass filter and 
becomes worse in dewhitening filter followed by prewhitening 
filter. This trend is the result of the relative magnitudes 
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of the tap outputs with respect to the noises caused by the 
device limitations. 

These problems encountered in the Reticon TAD-12 tapped 
delay line just showed that it is still in its developing 
stage. Improved results can be expected when the causes of 
the device limitations are better understood and corrected. 
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APPENDIX 



C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 



c 

c 

c 



c 

c 

c 



c 

c 



4 :^: ^ ^ 3 jc^ :$c:^:{c;jc:{c: 4 c 3 )(:!c : 4 i :4c :41:41:4c 

FFOGFAP OF C6SIGNING NQNRECURSIVE PREWHITENING FILTER BY 
VINCCW METHOD 

« :4t « :(i 3jc :4c :4c :4i :4! :jt 34 c :jt :4c * :4c 34 c :(t :4c =4t 4= 5»! ♦♦*♦♦* * 



THIS PROGRAM CALCULATES AND PLOTS AMPLITUDE ANC PHASE OF 
FIR CIGITAL FILTER WITH TRANSFER FUNCTION 
H(Z) = C(0)+C(1 )*Z*=!'(-l) + C(2)*Z=<'*(-2)+....+C(M )*2=i'*(-M) 
C(M) FINITE IMPULSE RESPONSE 

A SLOPE CF FREQUENCY AMPLITUDE RESPONSE FRCM FO TO 

HIGHEST FREQUENCY FP 

R CONSTANT FREQUENCY AMPLITUDE RSEPONSE FRCM 0 TO FG 

T = 1. (NORMALIZED DELAY UNIT) AND FP=0.5=<'( 1/T ) AND 
£ = FQ/FP 



COMPLEX ZINV,G,P 
COMPLEX CEXP.CMPLX 

DIMENSION F(150),C(30) ,W(30 ),D(30)tDD(30),CW(20), 
lAMPdSC )tPHASE( 150) 

INTEGER’»4 ITB ( 12 )/ 12*0/ 

REAL*A RTB(28)/28*0.0/ 

CATA R/0.1/,A/1./ ,E/119.048/,N/5/ 

PI=2. 14159265 

CCNVRT=180./PI 

M=2*N+1 

WFITE(6,500) M 

5QC FORMAT! IX, 'LENGTH OF FILTER IS:', 14///) 

TO COMPUTE FINITE IMPULSE RESPONSE ANC WINDOW FUNCTION 



W (N*l)=l. 

0 (N+1 )=R+A/(2*PI )*( .5-l./(2.*E) )**2 
CO 600 K=1.N 

C (K)=A* (CCS(K*PI)-CCS (K*PI/E) )/ ( 2.*PI**3*K = 



C (N + 



C (N + l+K ) = C(N+1-K) 



DC (K 



-K )=D(K) 



*2 ) 



=0.5^+0.46*CQS(2.*PI*K/M) 



W(N + 1-K )=DD(K) 

W (N + l+K ) = W(N+1-K) 

600 CONTINUE 

WRITE (6,650) (C( I) ,I=1,M) 

650 FCFMAT(//2X, 'COEFFICIENTS A RE : ' /8 ( IX , FIO. 5 ) ) 
CMAX=A6S(C(N+1 ) ) 

DC 700 1=1, M 

IF(ABS(C(I )).GT.CMAX) CMAX=A8S(C(D) 

700 CONTINUE 

CO 750 1 = 1, M 

NORMALIZE FIR 

C (I) = C( D/CMAX 
750 CONTINUE 

WFITE(6,800)(C( I) ,I=1,M) 

800 FCPMAT( //2X, 'NORMALIZED COEFFICIENTS AR E : • / £ ( IX, F 10. 5 ) 1 
WRITE (6,900) (W(I),I=1,M) 

900 FCFMAT(//2X, 'HAMMING WINDOW ARE : « /8( IX, FIO. 5) ) 

DC 9 1=1, M 
9 CW(I)=C (I )*W(I ) 

WFITE(6,1000) (CW( I ) ,I = 1,M ) 

1000 FORMAT! //2X,' COEFFICIENTS IMPROVED BY HAMMING ARE' 
l/8( 1X,F10.5) ) 

WFITE(6,50) 

50 FORMAT! ///5X,' FREQUENCY ',8X,* AMPLITUDE' , 13X , • PHAS 5 ' /// ) 
CO 10 J = 1,1Q1 



F IS FREQUENCY IN FRACTIONS OF NYQUIST FREQUENCY 



150 



c 

FU)=0,C05«PL0AT( J-1) 
2INV=CEXP(CMPLX(0.,-2.*FC J)*PI ) ) 

P = C^.PLX(1.,0.) 

G = Cf^PLX(0.,0.) 

CC 6 1 = 1, M 
G=G4p*CW(I ) 

6 P=F*Z!NV 

Af'P (J) = CABS(G) 

X=PEALCG) 

V=AIMAG(G» 

IF(ABS(X).GT.l.E-8) GO TO 60 
FF4SE(J)=SIGN(90.,Y> 

GC TO 8 

60 PFASE(J)=ATAN2(Y,X)*C0NVRT 
8 UPITE16,20) F( J) , ANP(J» ,PHASE( J ) 

20 FCPMAT(5X,F6.3,4X,E16.8,4X,E16.e> 

10 CCMINUE 

ECUVALENCE(TITLE,RTB(5) ) 

RE^L*8 TITLE(12) 

PEAC(5,300) TITLE 
300 FCPMAT(6A8) 

C^LL CRAWP(101,F,AMP,IT8,RTE) 

PE4C(5,300) TITLE 

C^LL DRAWPdOl, F, PHASE, ITB, RTB) 

STOP 

ENC 

//GC.SYSIN CC * 

AMPLITUCE FRECUENCY RESPONSE OF PREWHITENING 
FILTER / LENGTH = 11 / HARMING WINDOW 
PHASE RESPONSE OF PREWHITENING FILTER 
LENGTH * 11 / HAMMING WINDOW 
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APPENDIX E 



C FFGGFAM OF CESIGNING DWF FROM GIVEN IMPULSE RESPONSE OF 
C FWF EY FRECUENCY SAMPLING METHOD 

C TFE FPESCRIBEC CHARACTERISTIC AND DESIGNED FRECUENCY 
C CFAP/CTEPISTIC OF DWF ARE PLOTTED ON THE SAME GRAPH 

C 4= ’S' ^ 4: =!'=!' ’S' T =S= =4= * ’i' * ’i' # 



THIS PROGRAM CALCULATES AND PLOTS AMPLITUDE ANC PHASE OF 
FIR DIGITAL FILTER WITH TRANSFER FUNCTION 
H(2» = C(0 )+C(l)*Z**(-l ) + C(2 )*Z**(-2)+....+C(M )«Z**( -M ) 
C(M) FINITE IMPULSE RESPONSE 



CCHFLEX ZINV,G,P»H(150) 

COMPLEX CEXP.CMPLX 

DIMENSION F(150),C(30),W(30),D(20),0D(30)»CW(30), 
lAMP (1 50 )» PHASE ( 150 )»HH( 150) 

INTEGER=»4 ITB( 12 ) / 12*0/ 

REAL*4 RTB(28 )/28*0.0/ 

FACT0R=1.12 
MM=0 
N = 5 

PI=3.1415S265 
CCNVRT= 180-/PI 
M=2*N+1 

WFITE<6,500) M 

500 FCFMAT( IX, 'LENGTH OF FILTER IS:',IA///) 

READ(5, 1) (C(I) ,I=1,M) 

1 FCRMAT( EF1C.5) 

610 WRITE (6,650MC( I ) ,I=1,M ) 

650 FCFMAT( //2X, 'COEFFICIENTS ARE : ' /8 ( IX , FIO .5 ) ) 

CCMPLTE HAMMING WINDOW 

W (N+l)=l. 

DC 6C0 H=1,N 

CC(K)=0.54+0.46*C0S (2.*PI*K/M) 

WCN + l-K ) = 0D(K ) 

W (N + l+K )=W (N+l-K ) 

60C CONTINUE 

CMAX=ABS(C{N+1) ) 

DC 700 1=1, M 

IF(ABS( C(I) ).GT.CMAX) CMAX= ABS ( C ( I ) ) 

700 CONTINUE 

DC 750 1 = 1 ,M 

NORMALIZE FIR 

0 (I)=C( i)/CMAX 
750 CONTINUE 

WRITE (6,800MC( I ) ,I = 1,M) 

800 FCRMAT( //2X, 'NORMALIZED COEFFICIENTS ARE : * /8( IX, F10.5) J 
DC 5 1=1, M 

< cw(n=c (I )*w(i ) 

WPITE(6,1000)(CW( I ) ,I=1,M) 

1000 FCRMAT( //2X,'CCErFICIENTS IMPROVED BY HAMMING ARE' 
l/e( 1X,F10.5) ) 

WFITE(6,50) 

50 FC F MAT (///5X,'FREQUENCY',8X,' AMPLITUDE' , 13X , • PHAS E ' / // ) 
CO 10 J=l,101 

F IS FREQUENCY IN FRACTIONS OF NYQUIST FREQUENCY 

F( J )=0.Q05*FL0AT( J-1) 

ZINV=CEXP(CMPLX(0.,-2.*F( J )*PI ) ) 

F=CMPLX(1.,0.) 

G=CMFLX (0.,0.) 
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11 

e 



cc 6 1=1, M 

IF(MM.£C.l ) GO TO 11 
G=G+P*CW(I ) 

GC TO 6 

G = G4P*CW(n*FACT0R 

F=F*ZINV 

^^'F (J)=CABS(G) 



FPECteNCY SAMPLE POINTS OF DWF 

MJ) = 1./G 
HK J)=C4BS(H(J) ) 

X=REAL(G)- 

Y=AIMAG{G) 

IF(AeS(>).GT.l.E-8) GO TO 60 
PF^SE(J)=SIGN(90.,Y) 

GC TO 8 

6C PF^SE( J)=ATAN2(Y,X)*C0NVRT 
8 WPITE<6,20) F( J) ,AMP( J) ,PHAS£( J) 
2C FCFMAT(5X,F6.3,4X,£16.8,4X, E16.8) 
10 CONTINUE 

ECU IVALENCE (TITLE, RTB(5) ) 

RE4L*8 TITLE(12) 

IF(MM.EC.l) GO TO 301 
GC TO 302 
301 ITE(1)=2 
ITE(2)=0 

C4LL CRAWPdOl, F,AMP, ITB,RTB) 

202 L=N+1 



CF CWF 



1 = 1, L 



HOC 



1500 



300 



CC 1500 
H = I-1 

6G»CMPLX(0.,0. ) 

CC 1100 J=l,101 

2=CEXP(CMPLX(0.,2.*PI*F(J )#K) ) 
GG=6G+2#H( J) 

CCNTINUE 
C (n = 66/101. 

C (M + l-I ) = C(I ) 

CCNTINUE 



NN=MM+1 
IF(MM.EC.2» 
PE4C(5, 200) 
FCFMAT(6A8) 
ITE(1)=1 
ITE(2)=1 



GO TO 
TITLE 



2100 



FLCT THE PRESCRIBED FREQUENCY CHARACTERISTIC CF CWF 

LSING INVERSE DISCRETE FCURISR TRANSFORM TO CALCULATE THE 
IMPULSE RESPONSE OF CWF 

CALL CRAWP(101,F,HH,ITB,RTB ) 

GC BACK TO PLOT THE AMPLITUDE FREQUENCY RESPONSE OF DWF 

610 



GC TO 
2100 STCF 
ENC 

//GC.SYSIN CD 

-0.00091 0.0000 

-0.21^27 C.O 
AMPLITUDE FRECUENCY 
FILTER / LENGTH = 11 



-0.01219 
-0.01219 
RESPONSE OF 
/ HAMMING 



0.0000 
0.0 
DEWHITEMNG 
WINDOW 



-0.21A27 

-0.00091 



1.0000 
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APPENDIX C 



Q # :»5!!4e :je :{:5is * 4: Sr tS * 5js ’Jt ^ 5> ♦ Jf! ^ 4 ^ ^ * 

C* FFCGFAM OF DESIGNING NONFECURSIVE BANDPASS FILTER 

« jji *♦ s» ♦:» sjt » sfs ♦««« 4c !jt Jjt ## J)I :{e !}C Jji -t 



7HS PROGRAM CALCULATES AND PLOTS AMPLITUDE AND PHASE OF 

FIR DIGITAL FILTER WITH TRANSFER FUNCTION 

M2) = C(0)+C(l)*Z**(-l) + C(2)*Z**(-2)+....+C(M )«Z**(-M) 

C(M) FINITE IMPULSE RESPONSE 

A AMPLITUDE OF BANDPASS 

FI e F2 I CRITICAL FREQUENCIES 



CCMFLEX ZINVfGfP 
CCMFLEX CEXPrCMPLX 

DIMENSION F(150) »C(30) ,W( 30),D< 20) ,DD(30) tCM3C) , 
lAMPdSO ),PHASE(150) 

INTEGER*^ ITB< 12)/12’!'0/ 

REAL*4 RTB(28)/28*0.0/ 

DATA A/l./,Fl/0.1980/,F2/0.2020/,N/5/ 

PI=2. 14159265 

CCNVRT=180./PI 

M=2*N+1 

WP1TE(6,500) M 

5QC FCRMATdX, 'LENGTH OF FILTER IS:*, 14///) 



W (N+l) = l. 

C (N + 1)=2.*A«(F2-F1) 

DC 6C0 K=1,N 

D (K ) = A* (SIN(2.’*‘PI*K’!'F2)-SIN (2.*PI*K*FI) )/(K’»PI ) 

C (N + l-K )=D(K) 

C (N+1+K) = C(N+1-K) 

DC (K) =0.54+0. 46*C0S (2.*PI=>K/M) 

W(N + 1-K )=CC(K) 

W (N + l+K ) = W(N + 1-K) 

600 CONTINUE 

WRITE (6,650) (C( I ) ,I = l,M) 

650 FCRMAT(//2X, 'COEFFICIENTS ARE : ' / 8 ( IX, F 10. 5 ) ) 
CMAX=ABS(C(N+1) ) 

DC TOO 1=1, M 

IF(ABS(C (I ) ).GT.CMAX) CMAX=ABS ( C ( I ) ) 

700 CONTINUE 

DC 750 1=1, M 

NORMALIZE FIR 

C (I ) = C(I )/CMAX 
75C CONTINUE 

WRITE (6,800) (C( I) ,I=1,M) 

8QC FCRMAT(//'2X, 'NORMALIZED COEFFICIENTS AR E : ' / E ( IX, F 10 . 5 ) ) 
WRITE (6,900) (W(I),I=1,M) 

900 FORMAT (//2X, 'HAMMING WINDOW ARE : ' /8 ( IX , F 10. 5 ) ) 

DC 9 1=1, M 
9 CW(I)=C(I )*W(I ) 

WPITE(6,1000)(CW( I ) ,I=1,M) 

1000 FCRMAT(//2X, 'COEFFICIENTS IMPROVED BY HAMMING ARE* 
1/6(1X,F10.5)) 

WRITE(6,50) 

50 FORMAT ( ///5X , • F RE QUENCY * , 8X , • AMPL ITUDE ' , 13X , • PHAS E * /// 1 
DC 10 J = l, 101 

F IS FREQUENCY IN FRACTIONS OF NYQUIST FREQUENCY 

F(J)=0.01=»FL0AT(J-1) 

ZINV=CE)(P(CMPLX(0.,-2.=*‘F( J)*PI ) ) 

F=CMPLX(1.,0.) 

G=CMPLX (0.,0.) 
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CC 6 1=1, M 
G=G+P*CW(I ) 

F = F=»ZINV 
^NP(J)=CABS(G) 

X=FEAL(G) 

Y=AIMAG<G) 

IF(ABS(X).GT.l.E-8) GO TO 60 
PI-/SE{J) = SIGN(90.,Y) 

GC TO 8 

PFASE( J )=ATAN2( Y,X)«CONVRT 
V»PITE(6,20) F(J),A^'P(J),PHASE(J) 
FOPFAK 5X,F6.3,4X,E16.8,4X,E16.8 ) 
CCNTINUE 

ECUIVALENCE(TITLE,RTB(5) ) 

REAL*8 TITLE(12) 

REAC(5,20C) TITLE 
FCF^AT(6A8) 

CALL DRAWP(101,F,AMP,ITB,RTE) 
PEAC(5,30C) TITLE 
CALL DRAWP(10l,F, PHASE, ITBfPTB) 
STOP 
E^C 

//GC.SYSIN DC * 

AMPLITUDE FREQUENCY RESPONSE OF BANDPASS 
FILTER / LENGTH = 11 / HAMMING WINDOW 
PHASE RESFCNSE OF BANDPASS FILTER 
LENGTH = 11 / HAMMING WINDOW 



6C 

E 

20 

10 



iOQ 
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APPENDIX D 



C THIS PROGRAM CALCULATES IMPULSE RESPONSE GIVEN AMPLITLDE 
C FPEQLENCY RESPONSE OF FILTER 
COMPLEX CEXP.CMPLX 

CIMENSICN F (150) , COO), WOO ),CW (30 ,00(301 ,F( 150) 

N = 5 

FI=3.1A15S265 

F=2*N+1 

WFITE(6,5C0) M 

5CC FCPMAT( IX, 'LENGTH OF FILTER IS: *,IA//) 

COMPUTE HAMMING WINDOW 

W(N-H)=1. 

DC EOO K=1,N 

CC (K)=0.5 4+0.46«C0S (2. =<'PI *K/M ) 

W(N+1-K)=00(K) 

W(N+1+K ) = W(N-H-K) 

600 CONTINUE 

PEAC(5,'?00) (H( J) ,J = 1,101) 

700 FGFMAT(6F10*5) 

WRITE(6,750) ( H( J ) , J* 1, 10 1 ) 

75C F0RMAT(//2X,'DATA ARE :',///( 8 ( 1 X , F10.5 ),/) ) 

NCFMALIZE MJ) 

FMAX=H( 69) 

CC 2000 J=l,101 

IF(H( J) .GT.HMAX) HMAX=H(J) 

2000 CONTINUE 

CC 2100 J=l,101 
H(J)=H( Jl/HMAX 
2100 CONTINUE 

LSIN6 INVERSE DISCRETE FOURIER TRANSFORM 
TC CALCULATE THE IMPULSE RESPONSE 



L=M 

CC 1500 1=1, L 
K = I-1 

GC-=CMPLX(0,,0. ) 

DC 1100 J = l,101 
F(J)=0.01*FLOAT(J-1) 

Z=CEXP(CMPLX(0.,2.*PI«F(J)*K) ) 

GG = GG + Z=»H(J) 

HOC CONTINUE 

C (I ) = GG/1C1. 

1500 CONTINUE 

WPITE(6,650)(C(I) ,I=1,M) 

650 FCRMAT(//2X, 'COEFFICIENTS ARE : ' /8 ( IX , F 10. 5 ) ) 

NORMALIZE FIR- 

CMAX=ABS(C(N+1 ) ) 

CC 800 1=1, M 

IF(ABS( C(I n.GT.CMAX) CMAX=ABS ( C ( I ) ) 

800 CONTINUE 

CC 900 I=1,M 
C(I ) = C( D/CMAX 
900 CONTINUE 

WRITE (6,850) (C( I ) ,1 = 1 ,M) 

850 FCFMAT(//2X, 'NORMALIZED COEFFICIENTS AR E : • / £ ( IX , F 10. 5 ) 
CC 9 1=1, M 
9 CW(I)=C(I)=f=W(I) 

WRITE (6, 1000 ) (CW( I ) ,I = 1,M ) 

1000 FCRMAT( //2X, 'COEFFICIENTS IMPRCVEC BY HAMMING' 
l/8( 1X,F10.5) ) 

S7CP 

ENC 
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APPENDIX E 



C PPCGPAM FOR PLOTTING 100 CALCULATED AND MEASURED POINTS 
C CF ANY NONPECURSIVE FILTER WITH GIVEN IMPULSE RESPONSE 

C :Ji ♦ 41 4 * :Jr jjt :ji 4t ♦ :^c sjc 5jc :jt:Jc * sjt * * * * 5jt ^ :jc :jt :jc :{; 5jt * 4e :je :» >!t ;|t « 4e * * * ♦ * :{c :(! ♦ :(! 5{c 



This PROGRAM CALCULATES AND PLOTS AMPLITUDE AND PHASE QF 
FIR DIGITAL FILTER WITH TRANSFER FUNCTION 
M2) = C(0 )+C(l )*Z**(-l) + C(2)*Z**(-2)+....+C(M )*Z*=«'(-M) 
C(M) FINITE IMPULSE RESPONSE 



CCNPLEX ZINV,G,P 
CCNFLEX CEXP.CMPLX 

DIN ENSIGN F(150) ,C(30) ,AMP( 150) , PHASE (150) ,B(150) , 
166(150) iB A (150), AMPL (150), AMPLT( 150) 

1AVPL( 150) ,AMPLT(150) 

INTEGER^4 ITB ( 12 ) /12*0/ 

REAL*4 RTB(28 )/28*0.0/ 

PI=2. 1^159265 
N = ll 

CCNVRT=180./PI 
REAC(5, 500)(C( I) ,I=1,M) 

50C FCPNAT(8F10.5) 

REAC(5,500) (B(K) ,K=1,101) 

NPITE(6,850) (B(K) ,K=1, 101) 

85C FCPMAT(//2X,' DATA ARE : • ,/// (8 ( IX , F10.5 ) , / ) ) 

BMAX=B( E9) 

CC 1000 J=l,101 

1000 IF(B( J).GT.BMAX) BMAX=B(J) 

DC 1100 J=l,101 
EE(J)=B(J)/BMAX 
1100 6A( J)=20.*ALOG10(B6(J ) ) 

WRITE (6,80C)(C( I),I=1,M) 

800 FCPMAT(//2X, 'NORMALIZED COEFFICIENTS AR E : • / €( IX, F 10. 5 ) ) 
WRITE(6,50) 

50 FC P MAT (///5X, 'FREQUENCY', 8X , • AMPL ITUDE ’ , 13X ,' PHASE • /// ) 
CC 10 J=l,101 

F IS FREQUENCY IN FRACTIONS OF NYQUIST FREQUENCY 

F( J)=0.005*FLOAT( J-1) 

ZINV=CEXP(CMPLX(0.,-2.=!‘F(J )*PI ) ) 

F=CNPLX(1.,0.) 

G=CNPLX(0.,0.) 

CC 6 1=1, M 
G = C-4P*C (I) 

6 P = F’»ZINV 

AMP (J)=CABS(G) 

X=PEAL(G) 

Y=AIMAG(G) 

IF (ABS(X).GT.l.E-8) GO TO 60 
PHASE(J)=SIGN(90.,Y) 

GC TO 8 

60 PHASE! J)=ATAN2(Y,X )*CCNVRT 
6 WPITE(6,20) F( J) ,AMP(J) , PHASE! J ) ,E(J) 

20 FCFMAT(5X,F6.3,4X,E16.8,4X,E16.8,6X,E16. 8) 

10 CONTINUE 

AHAX=AMP( 51) 

CC 2C00 J=l,101 

IF (AMP( J ) .GT.AMAX ) AMAX = AMP(J) 

2000 CONTINUE 

CC 2100 J=l,101 
AMPL! J)=AMP( J)/AMAX 
2100 AMPLT! J )=20.*AL0G10(AMPL( J) ) 

IT8 a )=1 
ITE(2)=0 

ECUIVALENCE(TITLE,RTB(5) ) 
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RE4L*8 TITLEdZ) 

RE4C(5,7) TITLE 
7 FCFN4T(648) 

C4LL CR4WP(101,F,AMPLTtITB,RTB) 

ITE(1)=3 

ITE (2)=3 

CALL CRAWP(101,F,BA,TTB»RTB) 

STCP 

ENC 
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APPENDIX F 



TAPPING RESISTANCES FOR NON-RECURSIVE FILTER 




out 



PIGDRE P.1 Coefficient Implementation 



The contribution of a single tap,Vj^, to the output, 
is the product of and the coefficient, Aj^, or 




I I _ ^f 

I Vk I ■ r + Rk 



The expression indicates that Rk is inversely proportional 

to Ak> therefore the largest coefficient will correspond to 

the minimum resistance, R„-„, so 

min 

A = 

max r + R^in 



Solving this equation for R£ and substituting that result 
into the Aj^ expression yields: 

I ^ I _ \iax ^ ^in^ 

'^k' r + Rj^ 
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or, 

~ * Rmin) ■ r 

For the TAD- 12, the output impedance for a tap, r, is about 
5 Kilo-ohms. 
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